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Abstract. This set of lectures describes the physics of moment formation, the basic physics of the 
Kondo effect and the development of a coherent heavy electron fluid in the dense Kondo lattice. The 
last lecture discusses the open problem of quantum criticality in heavy electron systems. 
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1. LOCAL MOMENT FORMATION 
1.1. Introduction 

The last two decades have seen a growth of interest in "strongly correlated electron 
systems": materials where the electron interaction energies dominate the electron kinetic 
energies, becoming so large that they qualitatively transform the physics of the medium. 

m 

Examples of strongly correlated systems include 

• Cuprate superconductors, where interactions amongst electrons in localized 3d- 
shells form an antiferromagnetic Mott insulator, which develops high temperature 
superconductivity when doped. 

• Heavy electron compounds, where localized magnetic moments formed by rare 
earth or actinide ions transform the metal in which they are immersed, generating 
quasiparticles with masses in excess of 1000 bare electron masses.[|l]| 

• Fractional Quantum Hall systems, where the interactions between electrons in the 
lowest Landau level of a semi-conductor heteroj unction generate a new electron 
fluid, described by the Laughlin ground- state, with quantized fractional Hall con- 
stant and quasiparticles with fractional charge and statistics. ^ 

• "Quantum Dots", which are tiny pools of electrons in semiconductors that act as 
artificial atoms. As the gate voltage is changed, the Coulomb repulsion between 
electrons in the dot leads to the so-called "Coulomb Blockade", whereby electrons 
can be added one by one to the quantum dot. 

Strongly interacting materials develop "emergent" properties: properties which re- 
quire a new language[0 and new intellectual building blocks for their understanding. 
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FIGURE 1. Depicting localized 4/, 5/ and lid atomic wavefunctions. 



This chapter will illustrate and discuss one area of strongly correlated electron physics 
in which localized magnetic moments form the basic driving force of strong correlation. 
When electrons localize, they can form objects whose low energy excitations involve 
spin degrees of moment. In the simplest case, such "localized magnetic moments" are 
represented by a single, neutral spin operator 




where a denotes the Pauli matrices of the localized electron. Localized moments de- 
velop within highly localized atomic wavefunctions. The most severely localized wave- 
functions in nature occur inside the partially filled 4/ shell of rare earth compounds 
(Fig. P such as cerium (Ce) or Ytterbium {Yb). Local moment formation also occurs 
in the localized 5/ levels of actinide atoms as uranium and the slightly more delocal- 
ized 3d levels of first row transition metals(Fig. [1]). Localized moments are the origin 
of magnetism in insulators, and in metals their interaction with the mobile charge car- 
riers profoundly changes the nature of the metallic state via a mechanism known as the 
"Kondo effect". 

In the past decade, the physics of local moment formation has also reappeared in 
connection with quantum dots, where it gives rise to the Coulomb blockade phenomenon 
and the non-equilibrium Kondo effect. 



1.2. Anderson's Model of Local Moment Formation 



Though the concept of localized moments was employed in the earliest applications of 
quantum theory to condensed matter[], a theoretical understanding of the mechanism of 
moment formation did not develop until the early sixties, when experimentalists began 
to systematically study impurities in metals. Q 

In the early 1960s, Clogston, Mathias and collaborators!^ showed that when small 
concentrations n ■ of magnetic ions, such as iron are added to a metallic host, they develop 
a Curie component to the magnetic susceptibilty 

X = ^i^ M^=gj^lJ{J+\), (1) 

indicating the formation of a local moment. However, the local moment does not always 
develop, depending on the metallic host in which the magnetic ion was embedded. 
For example, iron dissolved at 1% concentration in pure Nb does not develop a local 
moment, but in the alloy Nb]^_^Mox a local moment develops for x > 0.4, rising to 
2.2/ig above x = 0.9. What is the underlying physics behind this phenomenon? 

Anderson[Q] was the first to identify interactions between localized electrons as the 
driving force for local moment formation. Earlier work by FriedelQ and Blandin[|T0|] 
had already identified part of the essential physics of local moments with the develop- 
ment of resonant bound-states. Anderson now included interactions to this picture. Much 
of the basic physics can be understood by considering an isolated atom with a localized 
S — 1/2 atomic state which we shall refer to as a localized "d-state". In isolation, the 
atomic bound state is stable and can be modeled in terms of a single level of energy 
and a Coulomb interaction 

U = ^J JWy(3c-/)|v/(x)|V(^)l^ (2) 

where V{x — j^) = e^/4-n£Q\x — x'\ is the Coulomb potential. 
The Anderson model for a localized impurity atom is given by 



H = Hc + H„.^ + H, + H^ (3) 



where = E^J^a^do describes an isolated atomic d-state of energy E^ and occu- 
pancy in the "up" and "down" state. = Uh,^n,, is the inter-atomic interaction 



Landau and Neel invoked the notion of the localized moment in their 1932 papers on antiferromag- 
netism, and in 1933, Kramers used this idea again in his theory of magnetic superexchange. 
^ It was not until the sixties that materials physicist could control the concentration of magnetic impurities 
in the parts per million range required for the study of individual impurities. Such control of purity evolved 
during the 1 950s, with the development of new techniques needed for semiconductor physics, such as zone 
refining. 



between the up and down d-electrons. The term 



4ri k ka ka 

describes the dispersion of electrons in the conduction sea which surrounds the ion, 
where c^^^ creates an electron of momentum k, spin o and energy e^. When the ion 
is embedded within a metal, the energy of the d-state is degenerate with band-electron 
states, and the term 

describes the hybridization that then takes place with the conduction electron sea, where 
ct^ a describes the creation of a d-electron. The matrix element of the ionic potential 
between a plane wave conduction state and the d-orbital is 

V. = {ko\V\d'a') = J d\e-'^'-%„{r)yr{r)5^^,. (4) 

where yA(r) is the wavefunction of the localized orbital and V-o„(?") is the ionic poten- 
tial. This matrix element will have the same symmetry as the localized orbital- a matter of 
some importance for real d-states, or f-states 0. However, for the discussion that follows, 
the detailed k dependence of this object can essentially be ignored . 
Let us first focus on the atomic part of H, 

a 

The four states of this ion are 

E{£)=2Ej + U 
\d^) Eld^)=0 (6) 

\d^^) \d^ i) E{d^)=E^. 

To obtain a magnetic doublet as the ground-state, the excitation energies out of the 
doublet state must be greater than zero, i.e 

E{d^)-E{d^)=E^ + U >0 Ej + U/2> ~U/2 

E{d^)-E{d^) = -E^>0 U/2>E^ + U/2 (7) 

so that for 

U/2>\Ej + U/2\, 



^ A direct calculation shows that 

V{k) ^47tr' J r^drj,{kr)V{r)R^{r) {1=2) (5) 

is the overlap of the radial wavefunctions /?p(r) of the d-state and the / = 2 partial wave state of the 
conduction electron, with the ionic potential. 
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FIGURE 2. Phase diagram for Anderson Impurity Model in the Atomic Limit. 

the isolated ion has a doubly degenerate magnetic ground-state, as illustrated in 0. We 
see that provided the Coulomb interaction U is large enough compared with the level 
spacing, the ground-state of the ion becomes magnetic. The d-excitation spectrum of the 
ion will involve two sharp levels, one at energy E^, the other at energy E^ + U . 

Suppose this ion is embedded in a metal: the free electron continuum is then pulled 
downwards by the work function of the metal so that now the d-level energy is degen- 
erate with conduction electron energy levels. In this situation we expect the d-level to 
hybridize with the conduction electron states, broadening the sharp d-level into a reso- 
nance with a width A = A(£'^^), where A(£) is given by Fermi's Golden Rule. 

A(e) = ;r£|y(X)|25(e^-e) = TtN{e)V^{£) (8) 

k 

where N{e) = 5(e^ — e) is the electron density of states (per spin). In the discussion 

that follows, let us assume that over the energy width of the resonance, V{e) and 
N{e) ~ A^(0) are essentially constant. 

When this hybridization is small compared with U , we expect the ground- states of 
the ion to be essentially that of the atomic limit. For weak interaction strength U the 
hybridization with the conduction sea will produce a single d-resonance of width A 
centered around E^. In Anderson's model for moment formation, when U > Uc ^ ttA 
the single resonance splits into two, so that for large U ^ Uc, there are two d-resonances 
centered around E^ and E^ + U , as shown in Fig. |3|. To illustrate the calculations 
that lead to this conclusion, let us use a Feynman diagram approach. We shall treat 
Hj — H^-^ + Hfj as a perturbation to the non-interacting part of the Hamiltonian to be 
H = H^-\-H^. The Green's functions of the bare d-electron and conduction electron are 
then denoted by 

---"""-> G^{k,i(On) = [icOn - %]"^ 

► G« = [{m,)-E,]-' 

whilst the Feynman diagrams for the hybridization and the interaction terms are then 




FIGURE 3. Illustrating how the d-electron resonance splits to form a local moment. A)U < nA, single 
half- filled resonance. B) U > TiA, up and down components of the resonance are split by an energy U. 
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-U (9) 
Quite generally, the propagator for the d-electrons can be written 

Gdai^)^[(0-E,-L,^i(0)]-' (10) 

where E^^(a)) is the the self-energy of the d-electron with spin a. We delineate between 
"up" and "down", anticipating Anderson's broken symmetry description of a local 
moment as a resonance immersed in a self-consistently determined Weiss field. The 
density of states associated with the d-resonance is determined by the imaginary part of 
the d-Green function: 

Pdai(o) = yrnG,^ico-i5). (11) 



The Anderson model for local moment formation is equivalent to the Hartree approxi- 
mation to the d-electron self-energy, denoted by 




The first term in this expression derives from the hybridization of the d-electrons with the 
conduction sea. Notice that the d-state fluctuates into all k-states of the conduction sea, 
so that there is a sum over k inside Z^{co). The second term is the Hartree approximation 
to the interaction self energy. We can identify the fermion loop here as the occupancy of 
the —a d-state, so that 

G,^(ft)) = [ft) - {E, + Un,_^) - El (ft))]-i (12) 

so the Hartree approximation is equivalent to replacing E^^ = E^ + Un^ The 
hybridization part of the self energy is 



Ei(ft) + /5)=£ 



(0-e-.-i5 
k k 

Notice that since 1/(x=f/5) = /'(I/jc) ±/;r5(jc), it follows that /mE^(ft)±/5) ==FA(ft)), 
so the imaginary part of this quantity has a discontinuity along the real axis equal to the 
hybridization width. Using (|8|), you can verify that we can now rewrite this as 



Ei(ft) + /5) 



n co-e-id 
de A(e) 

7t (0-£-id 



(13) 



Typically A(e) will only vary substantially on energies of order the bandwidth, so that 
over the width of the resonance we can replace A(e) — A. Moreover, for a broad band 
of width D, the real part ofE(ft)) ~ ^A(0)ln[(ft)-D)/(ft)-FD)] is of order ft)/D and can 
be ignored, or absorbed into into a small renormalization of E^. This allows us to make 
the replacement 

Ei(ft)±/5) =T«Asgn5 (24) 

so that 

The density of states described by the Green-function is a Lorentzian centered around 
energy E^^: 



moreover, the occupancy of the d-state is given by the d-occupation at zero temperature 
is 

- icot-(^) ^''^ 

This equation defines Anderson's mean-field theory. [| It is convenient to introduce an 
occupancy = Y.a'^da '''^^ magnetization M = n^^ — n^^, so that n^^^ = \ {n^ + oM) 
(a = ±1). The mean-field equation for the occupancy and magnetization are then 



M = - £ cTcot-W^^ 1] (16) 



To find the critical size of the interaction strength where a local moment develops, set 
A/ ^ 0+ (replacing the second equation by its derivative w.r.t. M), which gives 

= ^'^"^ ( a 

which can be written parametrically as 

^. + y = A(c + |(l-nJ(l+c2)) 

Vc = %lS,{\^-c^) (18) 

where c = cot j • The critical curve described by these equations is shown in Fig. |^. 

From the mean-field equations, it is easily seen that for = 1, when the d- levels 
are half filled, the critical value Uc = ttA. This enables us to qualitatively understand 
the experimentally observed formation of local moments. When dilute magnetic ions 
are dissolved into a metallic host, the formation of a local moment is dependent on 
whether the ratio U / nA is larger than, or smaller than zero. When iron is dissolved in 
pure niobium, the failure of the moment to form reflects the higher density of states and 
larger value of A in this alloy. When iron is dissolved in molybdenum, the lower density 
of states causes U > Uc, and local moments form. 



The quantity d(, — cot ' — j-^—^j is actually the phase shift for scattering an electron off the d- 

resonance (see exercise), and the identity n^^ = ^ 5^ is a particular realization of the "Friedel sum rule", 

which relates the charge bound in an atomic potential to the number of nodes (— introduced into 

the scattering state wavefunction. 



7.2.7. The Coulomb Blockade 



A modem context for the physics of local moments is found within quantum dots. 
A quantum dot is a tiny electron pool in a doped semi-conductor, small enough so that 
the electron states inside the dot are quantized, loosely resembling the electronic states 
of an atom. Unlike a conventional atom, the separation of the electronic states is of the 
order of milli-electron volts, rather than volts. The overall position of the quantum dot 
energy levels can be changed by applying a gate voltage to the dot. It is then possible 
to pass a small current through the dot by placing it between two leads. The differential 
conductance G = dl/dV is directly proportional to the density of states p{o)) inside the 
dot G p(0). Experimentally, when G is measured as a function of gate voltage Vg, the 
differential conductance is observed to develop a periodic structure, with a period of a 
few milli-electron volts, [j^] 

This phenomenon is known as the "Coulomb blockade" and it results from precisely 
the same physics that is responsible for moment formation. A simple model for a 
quantum dot considers it as a sequence of single particle levels at energies £^ , interacting 
via a single Coulomb potential U, according to the model 

Hdot = - 'VsK, + jN{N - 1) (19) 

where n-^^ is the occupancy of the spin o state of the A level, N = Y.xa'^Xa '^^'^^^ 
number of electrons in the dot and Vg the gate voltage. This is a simple generalization 
of the single atom part of the Anderson model. Notice that the capacitance of the dot is 
C^e'^/U. 

Provided that U is far greater than the energy separation of the individual levels, 
U » — £^,, the energy difference between the n electron and n + \ electron state 
of the dot is given hy E{n + I) — E{n) = nil — eVg. As the gate voltage is raised, the 
quantum dot fills each level sequentially, as illustrated in Fig. ^ and when eVg = U, 
the n-th level becomes degenerate with the Fermi energy of each lead. At this point, 
electrons can pass coherently through the resonance giving rise to a sharp peak in the 
conductance. At maximum conductance, the transmission and reflection of electrons is 
unitary, and the conductance of the quantum dot will reach a substantial fraction of the 
quantum of conductance, e^/h per spin. A simple calculation of the zero-temperature 
conductance through a single non-interacting resonance coupled symmetrically to two 
leads gives 

where the factor of two derives from two spin channels. At finite temperatures, the 
resonance becomes broadened by thermal excitation effects, giving 

,2 
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FIGURE 4. Variation of zero bias conductance G — dlj dV with gate voltage in a quantum dot. Coulomb 
interactions mean that for each additional electron in the dot, the energy to add one electron increases by 
U. When the charge on the dot is integral, the Coulomb interaction blocks the addition of electrons and 
the conductance is suppressed. When the energy to add an electron is degenerate with the Fermi energy 
of the leads, unitary transmission occurs, and for symmetric leads, G — 2e^/h. 



where /(e) = l/(e^^ + 1) is the Fermi function. When interactions are included, we 
must sum over the n-levels, giving (See Fig. H. ) 



G{eVg,T) = — d£ 



2^2 
h 



dfje-eVg) 
de 



{nU 



The effect of field on these results is interesting. When the number of electrons in 
the dot is even, the quantum dot is in a singlet state. When the number of electrons is 
odd, the quantum dot forms a local moment. In a magnetic field, the energy of the odd- 
electron dot is reduced, whereas the energy of the even spin dot is unchanged, with the 
result that at low temperatures 



E{2n + l)-E{2n) = InU - ^j^B 
E{ln + l)-E{ln+\) = {2n + 1)U + jlj^B 



(21) 



so that the vokages of the odd and even numbered peaks in the conductance develop an 
alternating field dependence. 

It is remarkable that the physics of moment formation and the "Coulomb blockade" 
operate in both artificial mesoscopic devices and naturally occurring magnetic ions. 



1.3. Exercises 

1. By expanding a plane wave state in terms of spherical harmonics: 

l^r\k) = e-'-' = 4Kl^i^j,ikr)Y*^{k)Y,Jr) 

l,m 

show that the overlap between a state with wavefunction {x\\i/) = R{r)Yi^^{f) with a 
plane wave is given by V{k) = {k\V\\lf) = V{k)Yii^{k) where 

V{k) = 4^i~' J drr^V{r)R{r)ji{kr) (22) 

2. (i) Show that 5 = cot"' i^t) ^^'^ scattering phase shift for scattering off a resonant 

level at position E^. 

(ii) Show that the energy of states in the continuum is shifted by an amount — AeS (e) / %, 
where Ae is the separation of states in the continuum. 

(iii) Show that the increase in density of states is given hy dd/dE = p^{E). (See chap- 
ter 3.) 



3. Derive the formula (20) for the conductance of a single isolated resonance. 



2. THE KONDO EFFECT 

Although Anderson's mean-field theory provided a mechanism for moment formation, 
it raised many new questions. One of its inadequacies is that of the magnetic moment is 
regarded as a broken symmetry order parameter. Broken symmetry is possible when the 
object that breaks the symmmetry involves a macroscopic number of degrees of freedom, 
but here, we are dealing with a single spin. There will always be a certain quantum 
mechanical amplitude for the spin to flip between an up and down configuration. This 
tunneling rate T^^ defines a temperature scale 

called the Kondo temperature, which sets the dividing line between local moment behav- 
ior, where the spin is free, and the low temperature limit, where the spin becomes highly 
correlated with the surrounding electrons. Experimentally, this temperature marks the 
low temperature limit of a Curie susceptibility. The physics by which the local moment 
disappears or "quenches" at low temperatures is closely analagous to the physics of 
quark confinement and it is named the "Kondo effect" after the Japanese physicist Jun 
Kondo. n 



The Kondo effect has a wide range of manifestations in condensed matter physics: 
not only does it govern the quenching of magnetic moments inside a metal, but it also 
is responsible for the formation of heavy fermion metals, where the local moments 
transform into composite quasiparticles with masses sometimes in excess of a thousand 
bare electron masses. [jT2|] Recently, the Kondo effect has also been observed to take 
place in quantum dots that carry a local moment. (Typically quantum dots with an odd 
number of electrons). [|3|] 

In this section we will first derive the Kondo model from the Anderson model, and 
then discuss the properties of this model in the language of the renormalization group . 



2.1. Adiabaticity 

Let us discuss some of the properties of the Anderson model at low temperatures 
using the idea of adiabaticity. We suppose that the interaction between electrons in the 
Anderson model is increased continuously to values U >> A, whilst maintaining the 
occupancy of the d- state equal to unity = 1. The requirement that n^=\ ensures that 
the d-electron density of states is particle-hole symmetric, which implies that E^ = and 
E'(0) = 0. 

When U >> A, we expect that the d-electron spectral function p^j = ^ImG^((0 — i5) 
will contain two peaks at CO = ±U /2. Since the total spectral weight integrates to unity, 
/ d(Op{(o) = 1, we expect that the weight under each of these peaks is approximately 
1/2. Remarkably, as we shall now see, the spectral function at ft) = is unchanged by 
the process of increasing the interaction strength and remains equal to its non-interacting 
value ^ 

p^(ft) = 0) = - 

This means that the d-spectral function must contain a narrow peak, of vanishingly small 
spectral weight Z << 1, height ^ and hence width A* = ZA << A. This peak in the 
d-spectral function is associated with the Kondo effect, and is known as the Abrikosov- 
Suhl, or the "Kondo" resonance. Let us see how this comes about as a consequence of 
adiabaticity. For a single magnetic ion, we expect that the interactions between electrons 
can be increased continuously, without any risk of instabilities, so that the excitations of 
the strongly interacting case remain in one-to-one correspondence with the excitations 
of the non-interacting case U = 0, forming a "local Fermi liquid". 

In this local Fermi liquid, one can divide the d-electron self-energy into two 
components- the first derived from hybridization, the second derived from interactions: 

E(ft)-z5) = iA + Zjico-iS) 

Zj{(0-i5) = (1-Z"i)ft) + /Aft)^ (23) 

The "wavefunction" renormalization Z is less than unity. The quadratic energy depen- 
dence of Sy(o;) ~ ft)^ follows from the quadratic energy dependence of the phase space 
for producing particle-hole pairs. Using this result, the form of the d-electron propagator 
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FIGURE 5. "Kondo resonance" in the d-spectral function. At large U for the particle hole-symmetric 
case where n^ — \, the d-spectral function contains two peaks about O) ±t/ /2, both of weight approx- 
imately 1 /2. However, since the spectral function is constrained by the requirement that pj\(i> = 0) = j, 
the spectral function must preserve a narrow peak of fixed height, but vanishingly small weight Z < < 1 . 



for at low energies is 



G,{co-id) 



1 



(0-iA-Zi{(0) 
Z 

(0-iZA-iO{(0^) 



(24) 



This corresponds to a renormalized resonance of reduced weight Z < 1, renormalized 
width ZA. One of the remarkable results of this line of reasoning, is the discovery that 
d-spectral weight 

Pjiw ~ 0) = -ImG^;(£0-/5)|^^o = - 

is independent of the strength of U. This result, first discovered by Langreth[jnj] guar- 
antees a peak in the d-spectral function at low energies, no matter how large U becomes. 
Since we also expect a peak in the d-spectral function around CO ~ ±U /2, this line of 
reasoning suggests that the structure of the d-spectral function at large U, contains three 
peaks. 



2.2. Schrieffer- Wolff transformation 



If a local moment forms within an atom, the object left behind is a pure quantum top- 
a quantum mechanical object with purely spin degrees of freedom. [] 

These spin degrees of freedom do interact with the surrounding conduction sea. In 
particular virtual charge fluctuations, in which an electron briefly migrates off, or onto 
the ion lead to spin-exchange between the local moment and the conduction sea. This 
induces an antiferromagnetic interaction between the local moment and the conduction 
electrons. To see this consider the two possible spin exchange processes 



e^ + d\ 
e^ + d\ 



AEj ~ 



U+E, 



(25) 



The first process passes via a doubly occupied singlet d-state, so it can only take place 
if the incoming conduction electron and d-electron are in a mutual 5 = state. In the 
second process, in order that the conduction electron can hybridize with the d-state, it has 
to arrive and depart in a state with precisely the same d- orbital symmetry. This means 
that the intermediate state formed in the second process must be spatially symmetric, 
and must therefore be a spin-antisymmetric singlet 5 = state. From these arguments, 
we see that spin exchange only takes place in the singlet channel, lowering the energy 
of the singlet configurations by an amount of order 



7 



1 1 

■ + 



1 



1 



E, + U 



(26) 
(27) 



where V is the size of the hybridization matrix element near the Fermi surface. If we 



introduce the electron spin density operator 5(0) = L 



where is the 



number of sites in the lattice, then we expect that the effective interaction induced by the 
virtual charge fluctuations will have the form 



where 5^ is the spin of the localized moment. Notice that the sign of J is 
antiferromagnetic. This kind of heuristic argument was ventured in Anderson's pa- 
per on local moment formation in 1961. The antiferromagnetic sign in this interaction 
was quite unexpected, for it had been tacitly assumed by the community that exchange 



^ In the simplest version of the Anderson model, the local moment is a S — 1/2, but in more realistic 
atoms much large moments can be produced. For example, an electron in a Cerium Ce^+ ion atom lives 
in a 4/' state. Here spin-orbit coupling combines orbital and spin angular momentum into a total angular 
moment y=/ — 1/2 = 5/2. The Cerium ion that forms thus has a spin j — 5/2 with a spin degeneracy 
of 2y + 1 =6. In multi-electron atoms, the situation can become still more complex, involving Hund's 
coupling between atoms. 



forces would induce a ferromagnetic interaction between the conduction sea and local 
moments. This seemingly innocuous sign difference has deep consequences for the 
physics of local moments at low temperatures, as we shall see in the next section. 

Let us now carry out the transformation a little more carefully, using the method of 
canonical transformations introduced by Schrieffer and Wolff[[T4|, [TS] ]. The Schrieffer- 
Wolff transformation is very close to the idea of the renormalization group and will 
help set up our renormalization group discussion. When a local moment forms, the 
hybridization with the conduction sea induces virtual charge fluctuations. It is therefore 
useful to consider dividing the Hamiltonian into two terms 

where A is an expansion parameter. Here, 



atomic 





" 








is diagonal in the low energy {Hj) and the high energy (fi or (f* {H^) subspaces, 
whereas the hybridization term 



jo 
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provides the off-diagonal matrix elements between these two subspaces. The idea of the 
Schrieffer Wolff transformation is to carry out a canonical transformation that returns 
the Hamiltonian to block-diagonal form, as follows: 



U 
























(28) 



This is a "renormalized" Hamiltonian, and the block-diagonal part of this matrix H* — 
P^H'Pj^ in the low energy subspace provides an effective Hamiltonian for the low energy 
physics and low temperature thermodynamics. If we set U = e^, where S = — is anti- 
hermitian and expand S in a power series 

then expanding using the identity e^Be^"^ = B+ [A, 5] + j^[A, [A, 5]] . . . 

e^(//i-FAr)e-^ = //i+A (r + [5i,//i]) Q[5i, [Si,//]] -f [Si,r] + [52,//i]^ +... 
so that to leading order 



and to second order 



e\H, + Xy)e-' = //i + A2 Q [S, , Y] + [S^,H^]^ +... 



(29) 



Since [5j , is block-diagonal, we can satisfy (|2^ ) to second order by requiring ^2 = 0, 
so that to this order, the renormalized Hamiltonian has the form (setting A = 1) 



where 



1 



is an interaction term induced by virtual fluctuations into the high-energy manifold. 
Writing 



'0 
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and substituting into (^, we obtain V = —sH^ + H^s. Now since {Hj)^^^ = £^5^/^ and 
i^H)ab = ^a^ab diagonal, it follows that 



V. 



ab 



ab 
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From (30), we obtain 



int) ab 



ab 



E'a-E^ 



(30) 
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Some important points about this result 



We recognize this result as a simple generalization of second-order perturbation 
theory which encompasses off-diagonal matrix elements. 



• H-^^^ can also be written 



where T is given by 



H. 
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int 



T{Ea) + T{E,) 



T{E) 
TabiE) 
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(31) 



is the leading order expression for the scattering T-matrix induced by scattering off 
'Y . We can thus relate H^^^^ to a scattering amplitude, and schematically represent it 
by a Feynman diagram, illustrated in Fig. |^. 

If the separation of the low and high energy subspaces is large, then the energy 
denominators in the above expression will not depend on the initial and final states 
a and so that this expression can be simplified to the form 
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FIGURE 6. T-matrix representation of interaction induced by integrating out high-energy degrees of 
freedom 



where AB^ = — is the excitation energy in the high energy subspace labeled 
by A, and the projector = • 

If we apply this method to the Anderson model, we have two high-energy subspaces, 
with excitation energies A£'(J^ — > cfi) = —E^ and AE{d^ d^) = E^ + U, so that the 
renormalized interaction is 



d^^d°+e- 



E, + U 



-E, 



Using the identity + o^^^ ■ — ^^ad^tc "^^^ ^^^^ renormalized Hamilto- 

nian in the form 

tiint = E h,V<^\a^<^k'fi-^d+^' 

d^+e-^d^ d^^d°+e- 
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+ 



1 



where 
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/3' 



K = 1) 



(32) 



(33) 



where we have replaced = 1 in the low energy subspace. Apart from a constant, the 
second term 



1 



'k.k'a 



1 



+ ■ 



1 



E, + U E, 



is a residual potential scattering term off the local moment. This term vanishes for 
the particle-hole symmetric case E^ = —{E^ + U) and will be dropped, since it does 
not involve the internal dynamics of the local moment. Summarizing, the effect of the 
high-frequency valence fluctuations is to induce an antiferromagnetic coupling between 
the local spin density of the conduction electrons and the local moment: 



^ = E^^^fccT^to + E V^^ta^^;^'^ • ^d 
ka Lk' 



(34) 



This is the infamous "Kondo model". For many purposes, the k dependence of the cou- 
pling constant can be dropped. In this case, the Kondo interaction can be written 
^int ~ ^^^(0)av/(0) ■ S^, where V^«(0) = -^E^^a electron operator at the origin 

and V/'''(0)av/(0) is the spin density at the origin. In this simplified form, the Kondo 
model takes the deceptively simple form 



^ = E ^k^\a^ka + -^V^' (0)0 ViO) ■ 5, . (35) 
ka 



In other words, there is a simple point-interaction between the spin density of the metal 
at the origin and the local moment. Notice how all reference to the fermionic character 
of the d-electrons has gone, and in their place, is a 5 = 1/2 spin operator. The fermionic 
representation ( |33| ) of the spin operator proves to be very useful in the case where the 
Kondo effect takes place. 



2.3. Renormalization concept 



To make further progress, we need to make use of the concept of renormalization. 
In a general sense, physics occurs on several widely spaced energy scales in condensed 
matter systems. We would like to distill the essential effects of the high energy atomic 
physics at electron volt scales on the low energy physics at millivolt scales without 
getting caught up in the fine details. An essential tool for this task is the "renormalization 
group". [|1§ |T7|, ig, 

The concept of the renormalization group permits us to describe complex condensed 
matter systems using simple models that reproduce only the relevant low energy physics 
of the system. The idea here is that only certain gross features of the high energy physics 
are relevant to the low energy excitations. The continuous family of model Hamiltonians 
with the same low energy excitation spectrum constitute a "universality class" of models. 
(Fig. ^ Suppose we parameterize each model Hamiltonian H{D) by its cutoff energy 
scale, D, the energy of the largest excitations. The scaling procedure, involves rescaling 
the cutoff D ^ D' = D/b where b > I, integrating out the excitations E E [D',D] to 
obtain an effective Hamiltonian //^for the remaining low-energy degrees of freedom. 
The energy scales are then rescaled, to obtain a new H{D') = bH^ Generically, the 
Hamiltonian will have the block-diagonal form 
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(36) 



where Hj^ and H^j act on states in the low-energy and high-energy subspaces respectively, 
and V and provide the matrix elements between them. The high energy degrees 
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FIGURE 7. Scaling concept. Low energy model Hamiltonians are obtained from the detailed original 
model by integrating out the high energy degrees of freedom. At each stage, the physics described by the 
model spans a successively lower frequency window in the excitation spectrum. 



of freedom may be "integrated out" f\ by carrying out a canonical transformation and 
projecting out the low-energy component i?^ 

H{D) ^ UH{D)U^ = 

By rescaling 

H{D') = bH^ (38) 

one arrives at a new Hamiltonian describing the physics on the reduced scale. The 
transformation from H{D) to H{D') is referred to as a "renormalization group" (RG) 
transformation. This term was coined long ago, even though the transformation does 
not form a real group, since there is no inverse transformation. Repeated application of 
the RG procedure leads to a family of Hamiltonians H{D). By taking the limit b ^ I, 



H 



H 



(37) 



The term "integrating out" is originally derived from the path integral formulation of the renormalization 
group, in which high energy degrees of freedom are removed by integrating over these variables inside 
the path integral. 



these Hamiltonians evolve continuously with D. Typically, H will contain a series of 
dimensionless coupling constants {g-} which denote the strength of various interaction 
terms in the Hamiltonian. The evolution of these coupling constants with cut-off is given 
by a scaling equation, so that for the simplest case 

dlnD ^^^^^'^^ 

A negative /3 function denotes a "relevant" coupling constant which grows as the cut- 
off is reduced. A positive /3 function denotes an " irrelevant " coupling constant which 
diminishes as the cut-off is reduced. There are two types of event that can occur in such 
a scaling procedure (Fig. ||): 

• i) A crossover. When the cut-off energy scale D becomes smaller than the charac- 
teristic energy scale of a particular class of high frequency excitations, then at lower 
energies, these excitations may only occur via a virtual process. To accommodate 
this change, the Hamiltonian changes its structure, acquiring additional terms that 
simulate the effect of the high frequency virtual fluctuations on the low energy 
physics. The passage from the Anderson to the Kondo model is an example of one 
such cross-over. In the renormalization group treatment of the Anderson model, 
when the band- width of the conduction electrons becomes smaller than the energy 
to produce a valence fluctuation, a cross-over takes place in which real charge fluc- 
tuations are eliminated, and the physics at all lower energy scales is described by 
the Kondo model. 

• ii) Fixed Point. If the cut-off energy scale drops below the lowest energy scale 
in the problem, then there are no further changes to occur in the Hamiltonian, 
which will now remain invariant under the scaling procedure (so that the /3 function 
of all remaining parameters in the Hamiltonian must vanish). This "Fixed Point 
Hamiltonian" describes the essence of the low energy physics. 



2.4. "Poor Man" Scaling 

We shall now apply the scaling concept to the Kondo model. This was originally 
carried out by Anderson and Yuval using a method formulated in the time, rather than 
energy domain. The method presented here follows Anderson's " Poor Man's" scaling 
approach, in which the evolution of the coupling constant is followed as the band-width 
of the conduction sea is reduced. The Kondo model is written 

m £ c\^a^pc,p.l (39) 

where the density of conduction electron states p(e) is taken to be constant. The Poor 
Man's renormalization procedure follows the evolution of J{D) that results from reduc- 
ing D by progressively integrating out the electron states at the edge of the conduction 



H = 



band. In the Poor Man's procedure, the band-width is not rescaled to its original size 
after each renormalization, which avoids the need to renormalize the electron operators 
so that instead of Eq. @, H{D') = H^^. 

To carry out the renormalization procedure, we integrate out the high-energy spin 
fluctuations using the t-matrix formulation for the induced interaction H-^^, derived in 
the last section. Formally, the induced interaction is given by 



where 



Xe\H) 



aX Ah 

E-El 



where the energy of state |A) lies in the range \p' ,D\. There are two possible intermedi- 
ate states that can be produced by the action of H^^"^ on a one-electron state: (I) either the 
electron state is scattered directly, or (II) a virtual electron hole-pair is created in the in- 
termediate state. In process (I), the T-matrix can be represented by the Feynman diagram 




o a' a' 

for which the T-matrix for scattering into a high energy electron state is 
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E-D 



k" 



la c<b\ 



>a c<b\ 



(40) 



In process (II), 




the formation of a particle-hole pair involves a conduction electron line that crosses 
itself, leading to a negative sign. Notice how the spin operators of the conduction sea 
and antiferromagnet reverse their relative order in process II, so that the T-matrix for 
scattering into a high-energy hole-state is given by 



k'pa'-kacy 



e^„e[-D-D+8D] 



1 



-J^pdD 



E-D 



(41) 



where we have assumed that the energies and e^, are negligible compared with D. 
Adding (Eq. 0) and (Eq. gl|) gives 



J^p5D 



£) Pa a'a 



(42) 



In this way we see that the virtual emission of a high energy electron and hole generates 
an antiferromagnetic correction to the original Kondo coupling constant 

J{D')=J{D) + 2J^p^-^ 

High frequency spin fluctuations thus antiscreen the antiferrromagnetic interaction. If 
we introduce the coupling constant g = pj, we see that it satisfies 



dg 
dlnD 



I5{g) = -2g^ + 0ig'). 



This is an example of a negative /3 function: a signature of an interaction which is weak 
at high frequencies, but which grows as the energy scale is reduced. The local moment 
coupled to the conduction sea is said to be asymptotically free. The solution to this 
scaling equation is 



l-2goHD/D') 



and if we introduce the scale 



Dexp 
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^go 



(43) 



(44) 



we see that this can be written 



2g{D') 



HD'/T^) 



This is an example of a running coupling constant- a coupling constant whose strength 
depends on the scale at which it is measured. (See Fig. 

Were we to take this equation literally, we would say that g diverges at the scale 
D' — T^. This interpretation is too literal, because the above scaling equation has only 
been calculated to order g^, nevertheless, this result does show us that the Kondo 
interaction can only be treated perturbatively at energy scales large compared with the 
Kondo temperature. We also see that once we have written the coupling constant in terms 
of the Kondo temperature, all reference to the original cut-off energy scale vanishes from 
the expression. This cut-off independence of the problem is an indication that the physics 
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FIGURE 8. Schematic illustration of renormalization group flow from a repulsive "weak coupling 
fixed point, via a crossover to an attractive "strong coupling" fixed point. 



of the Kondo problem does not depend on the high energy details of the model: there is 
only one relevant energy scale, the Kondo temperature. 

It is possible to extend the above leading order renormalization calculation to higher 
order in g. To do this requires a more systematic method of calculating higher order 
scattering effects. One tool that is particularly useful in this respect, is to use the 
Abrikosov pseudo-fermion representation of the spin, writing 

= 1. (45) 

This has the advantage that the spin operator, which does not satisfy Wick's theorem, is 
now factorized in terms of conventional fermions. Unfortunately, the second constraint 
is required to enforce the condition that 5^ = 3/4. This constraint proves very awkward 
for the development of a Feynman diagram approach. One way around this problem, is 
to use the Popov trick, whereby the d-electron is associated with a complex chemical 
potential 

T 

The partition function of the Hamiltonian is written as an unconstrained trace over the 
conduction and pseudofermion Fock spaces. 



Z = Tr 



g-/3(//+,;rf(n,-l)) 



(46) 



Now since the Hamiltonian conserves n^, we can divide this trace up into contributions 
from the d*^, d^ and d^ subspaces, as follows: 



Z = e'^'/^Zid^) +Z{d^) + e-'''/^Z{d^) 

But since 5^ = in the d^and d^ subspaces, Z{d^) = Z{d^) so that the contributions to 
the partition function from these two unwanted subspaces exactly cancel. You can test 
this method by applying it to a free spin in a magnetic field, (see exercise) 

By calculating the higher order diagrams shown in fig ^ , it is straightforward, though 
laborious to show that the beta-function to order is given by 

^ = m = -2/ + 2/ + 0{g') (47) 



FIGURE 9. Diagrams contributing to the third-order term in the beta function. A "crossed" propagator 
hne indicates that the contribution from high-energy electrons with energies |e^| g[D— 5D,D] is taken 
from this hne. 



One can integrate this equation to obtain 

A better estimate of the temperature where the system scales to strong coupling is 
obtained by setting D' = 7^ and g = 1 in this equation, which gives 

where for convenience, we have absorbed a factor into the cut-off, writing D = 
D^f . Thus, 

Tj, = Dy/2foe~^o (49) 

up to a constant factor. The square-root pre-factor in is often dropped in qualitative 
discussion, but it is important for more quantitative comparison. 



2.5. Universality and the resistance minimum 

Provided the Kondo temperature is far smaller than the cut-off, then at low energies 
it is the only scale governing the physics of the Kondo effect. For this reason, we expect 
all physical quantities to be expressed in terms of universal functions involving the ratio 
of the temperature or field to the Kondo scale. For example, the the susceptibility 

X{T) = ^F{^), (50) 

and the resistance 

\{T) = ^nY) (51) 

both display universal behavior. 



We can confirm the existence of universality by examining these properties in the 
weak coupling limit, where T » T^. Here, we find 
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27tJ^pS{S+\X, (5=-) 
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where n- is the density of impurities. Scaling implies that at lower temperatures Jp 
Jp + 2(7p)^ln ^, so that to next leading order we expect 



l(r) 

X{T) 



n.jS{S+\)[Jp + 2{Jp)Hn^f, 
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\-2Jp-A{Jpf\nj + 0{{Jpf) 



(52) 
(53) 



results that are confirmed from second-order perturbation theory. The first result was 
obtained by Jun Kondo. Kondo was looking for a consequence of the antiferromagnetic 
interaction predicted by the Anderson model, so he computed the electron scattering 
rate to third order in the magnetic coupling. The logarithm which appears in the electron 
scattering rate means that as the temperature is lowered, the rate at which electrons 
scatter off magnetic impurities rises. It is this phenomenon that gives rise to the famous 
Kondo "resistance minimum" . 

Since we know the form of T^, we can use this result to deduce that the weak coupling 
limit of the scaling forms. If we take equation (p^, and replace the cut-off by the 
temperature D ^ T, and replace go by the running coupling constant go 8{T), we 
obtain 

siT) = (54) 

which we may iterate to obtain 

2^(7^) = 



21n(^)+ln2g(r) 



ln(4 
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21n I ^ 



(55) 



Using this expression to make the replacement Jp g{T) in ( p^ and (]53|), we obtain 
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(56) 



(57) 



From the second result, we see that the electron scattering rate has the scale-invariant 
form 

l(r) = ^^(r/r^). (58) 



where ^(jc) is a universal function. The pre-factor in the electron scattering rate is 
essentially the Fermi energy of the electron gas: it is the "unitary scattering" rate, 
the maximum possible scattering rate that is obtained when an electron experiences a 
resonant njl scattering phase shift. From this result, we see that at absolute zero, the 
electron scattering rate will rise to the value \{J') = (0), indicating that at strong 
coupling, the scattering rate is of the same order as the unitary scattering limit. We shall 
now see how this same result comes naturally out of a strong coupling analysis. 



2.6. Strong Coupling: Nozieres Fermi Liquid Picture of the Kondo 

Ground-state 

The weak-coupling analysis tells us that at scales of order the Kondo temperature, 
the Kondo coupling constant g scales to a value of order 0(1). Although perturbative 
renormalization group methods can not go past this point, Anderson and Yuval pointed 
out that it is not unreasonable to suppose that the Kondo coupling constant scales 
to a fixed point where it is large compared to the conduction electron band-width 
D. This assumption is the simplest possibility and if true, it means that the strong- 
coupling limit is an attractive fixed point, being stable under the renormalization group. 
Anderson and Yuval conjectured that the Kondo singlet would be paramagnetic, with 
a temperature independent magnetic susceptibility and a universal linear specific heat 
given by Cy = Jky' temperatures. 

The first controlled treatment of this cross-over regime was carried out by Wilson 
using a numerical renormalization group method. Wilson's numerical renormalization 
method was able to confirm the conjectured renormalization of the Kondo coupling 
constant to infinity. This limit is called the "strong coupling" limit of the Kondo problem. 
Wilson carried out an analysis of the strong-coupling limit, and was able to show that 
the specific heat would be a linear function of temperature, like a Fermi liquid. Wilson 
showed that the linear specific heat could be written in a universal form 

Cy = yr, 

0.4128 ±0.002 
^ = T 87j ^''^ 

Wilson also compared the ratio between the magnetic susceptibility and the linear 
specific heat with the corresponding value in a non-interacting system, computing 

^ _ X/X^ ^X( \ ^ 2 (60) 

within the accuracy of the numerical calculation. 

Remarkably, the second result of Wilson's can be re-derived using an exceptionally 
elegant set of arguments due to Nozieres [p^ that leads to an explicit form for the strong 
coupling fixed point Hamiltonian. Nozieres began by considering an electron in a one- 
dimensional chain as illustrated in Fig. [1^. The Hamiltonian for this situation is 
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FIGURE 10. Illustrating the strong-coupling limit of the Kondo model 



tiMtice = E [c\U+^)caU) +H.c] + J c\{Q)a ^^c ^{Q) ■ S (61) 

Nozieres argued that the strong coupling fixed point will be described by the situation 

J » t.ln this limit, the kinetic energy of the electrons in the band can be treated as a 
perturbation to the Kondo singlet. The local moment couples to an electron at the origin, 
forming a "Kondo singlet" denoted by 

\Gs) = ^(\n)-\m (62) 

where the thick arrow refers to the spin state of the local moment and the thin arrow 
refers to the spin state of the electron at site 0. Any electron which migrates from site 
1 to site will automatically break this singlet state, raising its energy by 3//4. This 
will have the effect of excluding electrons (or holes) from the origin. The fixed point 
Hamiltonian must then take the form 

^lattice = L [c^cj(7+l)ccj(7)+H.c]+ weak interaction (63) 

where the second-term refers to the weak-interactions induced in the conduction sea by 
virtual fluctuations onto site 0. If the wavefunction of electrons far from the impurity 
has the form \f/{x) ~ sm{kpx), where kp is the Fermi momentum, then the exclusion of 
electrons from site 1 has the effect of phase- shifting the electron wavefunctions by one 
the lattice spacing a, so that now \l^{x) ~ sm{kpX — 5) where 6 = kp-a. But if there is one 
electron per site, then 2{2kpa/ {27u)) = 1 by the Luttinger sum rule, so that kp = 7c/ (2a) 
and hence the Kondo singlet acts as a spinless, elastic scattering center with scattering 
phase shift 

5 = 7t/2. (64) 

The appearance of 5 = ;r/2 could also be deduced by appealing to the Friedel sum 
rule, which states that the number of bound-electrons at the magnetic impurity site is 



£cr-2|±i = lb 1 11, SO that 5 = ;r/2. By considering virtual fluctuations of electrons 
between site 1 and 0, Nozieres argued that the induced interaction at site 1 must take 
the form 

because fourth order hopping processes lower the energy of the singly occupied state, 
but they do not occur for the doubly occupied state. This is a repulsive interaction 
amongst the conduction electrons, and it is known to be a marginal operator under the 
renormalization group, leading to the conclusion that the effective Hamiltonian describes 
a weakly interacting "local" Fermi liquid. 

Nozieres formulated this local Fermi liquid in the language of an occupancy- 
dependent phase shift. Suppose the ka scattering state has occupancy n^^, then the the 
ground-state energy will be a functional of these occupancies ^[{n^^}]. The differential 
of this quantity with respect to occupancies defines a phase shift as follows 

55 Ac 

5;^^ = ^;^-— (66) 

The first term is just the energy of an unscattered conduction electron, while 
5({n^^,},ej^) is the scattering phase shift of the Fermi liquid. This phase shift can 
be expanded 

^ k 

where the term with coefficient 4> describes the interaction between opposite spin states 
of the Fermi liquid. Nozieres argued that when the chemical potential of the conduction 
sea is changed, the occupancy of the localized d state will not change, which implies 
that the phase shift is invariant under changes in /i. Now under a shift the change 
in the occupancy Lfe^'^^cr — >^ ^MP, so that changing the chemical potential modifies the 
phase shift by an amount 

A5 = (a + ^>p)A;U = (68) 

so that a = — p^>. We are now in a position to calculate the impurity contribution to 
the magnetic susceptibility and specific heat. First note that the density of quasiparticle 
states is given by 

1^5 a 
dE Tt de Tt 

so that the low temperature specific heat is given by Cy = (7^„;j^ + Jj) where 

r. = \-^)-^ (70) 

where the pref actor "2" is derived from the spin up and spin-down bands. Now in a 
magnetic field, the impurity magnetization is given by 

M=^-^ (71) 



Since the Fermi energies of the up and down quasiparticles are shifted to Ep^ — > — oB, 
we have = opB, so that the phase-shift at the Fermi surface in the up and down 

scattering channels becomes 

^ k 
= — + ac75-4>pc75 

= ^ + 2aa5 (72) 

so that the presence of the interaction term doubles the size of the change in the phase 
shift due to a magnetic field. The impurity magnetization then becomes 

M, = xfi = 2 (^^^ (73) 

where we have reinstated the magnetic moment of the electron. This is twice the value 
expected for a "rigid" resonance, and it means that the Wilson ratio is 

W^ = 4^ = 2 (74) 

2.7. Experimental observation of Kondo effect in real materials and 

quantum dots 

Experimentally, there is now a wealth of observations that confirm our understanding 
of the single impurity Kondo effect. Here is a brief itemization of some of the most 



important observations. (Fig. 11.) 



A resistance minimum appears when local moments develop in a material. For ex- 
ample, in Nb^_^Mox alloys, a local moment develops for x > 0.4, and the resistance 
is seen to develop a minimum beyond this point. [|20|, ^A\\ 

Universality seen in the specific heat Cy = yF{T /Tj^) of metals doped with dilute 
concentrations of impurities. Thus the specific heat of Cu — Fe (iron impurities 
in copper) can be superimposed on the specific heat of Cu — Cr, with a suitable 
rescaling of the temperature scale. 



Universality is observed in the differential conductance of quantum dots[|25|, |2q ] 
and spin-fluctuation resistivity of metals with a dilute concentration of 
impurities.[|24|] Actually, both properties are dependent on the same thermal 
average of the imaginary part of the scattering T-matrix 

p,- = n-— [ dco 2Im[T(co)] 
' ' m J \ dco J 

G = ^ [ dcoi -^]7tpIm[T{co)]. (75) 
n J \ oco 
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FIGURE 11. (a) Sketch of resistance minimum in MoxNb^ _^ (b) Sketch of excess resistivity associated 
with scattering from an impurity spin. Right hand-scale- differential conductivity of a quantum dot. 



Putting Kp jdG)(-^^^T{G)) = t{co/Tj^,T/Tj^), we see that both properties 
have the form 



p = n^^^t(T/Tj.) 

^' 'nmp ^ ' ^' 

G = -Yt{T/T^) 



(76) 



where t{T /T^) is a universal function. This result is bom out by experiment. 



2.8. Exercises 



1. Generalize the scaling equations to the anisotropic Kondo model with an anisotropic 
interaction 

H,= E J"^\a<^c,p-S\ (77) 

|etl.|ei/.a=(-<:.:>',z) 

and show that the scahng equations take the form 

dJa 



d\nD 



-2J,J,p + 0{J'), 



where and {a,b,c) are a cyclic permutation of {x,y,z). Show that in the special case where 
Jx = Jy = the scaling equations become 



dlnD 
d\nD 



-2J,J^p + 0{j'), 

-2{jfp + 0{j'), (78) 



so that -~J± = constant. Draw the corresponding scaling diagram. 

2. Consider the symmetric Anderson model, with a symmetric band-structure at half filling. 
In this model, the d'^ and d^ states are degenerate and there is the possibility of a "charged 
Kondo effect" when the interaction U is negative. Show that under the "particle-hole" 
transformation 

'ki - d^^-d\ (79) 

the positive U model is transformed to the negative U model. Show that the spin operators 
of the local moment are transformed into Nambu "isospin operators" which describe the 
charge and pair degrees of freedom of the d-state. Use this transformation to argue that 
when U is negative, a charged Kondo effect will occur at exactly half-filling involving 
quantum fluctuations between the degenerate d'^ and d^ configurations. 

3. What happens to the Schrieffer- Wolff transformation in the infinite U limit? Rederive 
the Schrieffer- Wolff transformation for an N-fold degenerate version of the infinite U 
Anderson model. This is actually valid for Ce and Yb ions. 

4. Rederive the Nozieres Fermi hquid picture for an SU (N) degenerate Kondo model. Explain 
why this picture is relevant for magnetic rare earth ions such as Ce^~^ or Yb^~^. 

5. Check the Popov trick works for a magnetic moment in an external field. Derive the 
partition function for a spin in a magnetic field using this method. 

6. Use the Popov trick to calculate the T-matrix diagrams for the leading Kondo renormahza- 
tion diagramatically. 



3. HEAVY FERMIONS 



Although the single impurity Kondo problem was essentially solved by the early sev- 
enties, it took a further decade before the physics community was ready to accept the 




FIGURE 12. Illustrating how the polarization of spin around a magnetic impurity gives rise to Friedel 
oscillations and induces an RKKY interaction between the spins 



notion that the same phenomenon could occur within a dense lattice environment. This 
resistance to change was rooted in a number of popular misconceptions about the spin 
physics and the Kondo effect. 

At the beginning of the seventies, it was well known that local magnetic moments 
severely suppress superconductivity, so that typically, a few percent is all that is required 
to destroy the superconductivity. Conventional superconductivity is largely immune 
to the effects of non-magnetic disorder [] but highly sensitive to magnetic impurities, 
which destroy the time-reversal symmetry necessary for s-wave pairing. The arrival 
of a new class of superconducting material containing dense arrays of local moments 
took the physics community completely by surprise. Indeed, the first observations of 
superconductivity in UBe^j,, made in 1973 [ ^7| ] were dismissed as an artifact and had 
to await a further ten years before they were revisited and acclaimed as heavy fermion 
superconductivity. [E^, E^] 



Normally, local moment systems develop antiferromagnetic order at low tempera- 
tures. When a magnetic moment is introduced into a metal it induces Friedel oscillations 
in the spin density around the magnetic ion, given by 

{M{x)) = -Jx{x-i){S{^)) 
where J is the strength of the Kondo coupling and 

Xiq) = (80) 

k k+q k 

is the the non-local susceptibility of the metal. If a second local moment is introduced 
at location x, then it couples to {M{x)) giving rise to a long-range magnetic interaction 



' Anderson argued in his "dirty superconductor theorem" that BCS superconductivity involves pairing 
of electrons in states that are the time-reverse transform of one another. Non-magnetic disorder does not 
break time reversal symmetry, and so the one particle eigenstates of a dirty system can still be grouped 
into time-reverse pairs from which s-wave pairs can be constructed. For this reason, s-wave pairing is 
largely unaffected by non-magnetic disorder. 



called the "RKKY"[|3g] interaction, Q 



Hrkky = -J^X - /) S{x) ■ S{x'). (81) 

The sharp discontinuity in the occupancies at the Fermi surface produces slowly decay- 
ing Friedel oscillations in the RKKY interaction given by 

jRKKAr)--j'p^^ (82) 

where p is the conduction electron density of states and r is the distance from the 
impurity, so the RKKY interaction oscillates in sign, depending on the distance between 
impurities. The approximate size of the RKKY interaction is given by E^i^^^y ~ -^^P- 

Normally, the oscillatory nature of this magnetic interaction favors the development 
of antiferromagnetism. In alloys containing a dilute concentration of magnetic transi- 
tion metal ions, the RKKY interaction gives rise to a frustrated, glassy magnetic state 
known as a spin glass in which the magnetic moments freeze into a fixed, but random 
orientation. In dense systems, the RKKY interaction typically gives rise to an ordered 
antiferromagnetic state with a Neel temperature ~ J^p . 

In 1976 Andres, Ott and Graebner discovered the heavy fermion metal CeAl^. pT| ] 
This metal has the following features: 

• A Curie susceptibility ~ T at high temperatures. 

• A paramagnetic spin susceptibility x ~ constant at low temperatures. 

• A linear specific heat capacity Cy = yT, where 7~ l600mJ / mol / is approxi- 
mately 1600 times larger than in a conventional metal. 

• A quadratic temperature dependence of the low temperature resistivity p = Po + 

Andres, Ott and Grabner pointed out that the low temperature properties are those of 
a Fermi liquid, but one in which the effective masses of the quasiparticles are approx- 
imately 1000 larger than the bare electron mass. The Fermi liquid expressions for the 
magnetic susceptibility x and the linear specific heat coefficient y are 



X = (M. 



7 = ^N*iO) (83) 

where A^* (0) = ^N{0) is the renormalized density of states and Fq is the spin-dependent 
part of the s-wave interaction between quasiparticles. What could be the origin of this 
huge mass renormalization? Like other Cerium heavy fermion materials, the Cerium 
atoms in this metal are in a Ce^^{4-f^) configuration, and because they are spin-orbit 



named after Ruderman, Kittel, Kasuya and Yosida 




FIGURE 13. (a) Single impurity Kondo effect builds a single fermionic level into the conduction sea, 
which gives rise to a resonance in the conduction electron density of states (b) Lattice Kondo effect builds 
a fermionic resonance into the conduction sea in each unit cell. The elastic scattering off this lattice of 
resonances leads to formation of a heavy electron band, of width 7^. 



coupled, they form huge local moments with a spin of 7 = 5/2. In their paper, Andres, 
Ott and Graebner suggested that a lattice version of the Kondo effect might be responsi- 
ble. 

This discovery prompted Sebastian Doniach[]^] to propose that the origin of these 
heavy electrons derived from a dense version of the Kondo effect. Doniach proposed that 
heavy electron systems should be modeled by the "Kondo-lattice Hamiltonian" where a 
dense array of local moments interact with the conduction sea. For a Kondo lattice with 
spin 1/2 local moments, the Kondo lattice Hamiltonian takes the form 



H 



k ka ka 



ka 



ka 



k'p 



(84) 



Doniach argued that there are two scales in the Kondo lattice, the Kondo temperature 
mdEj^j^j^y, given by 



'^RKKY 



(85) 



When Jp is small, then Ej^j^j^y » T^, and an antiferromagnetic state is formed, but 
when the Kondo temperature is larger than the RKKY interaction scale, Tj.» Ej^j^^^y, 
Doniach argued that a dense Kondo lattice ground-state is formed in which each site 
resonantly scatters electrons. Bloch's theorem then insures that the resonant elastic 
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Tji ~ Dexp[—1/Jp] 
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FIGURE 14. Doniach diagram, illustrating the antiferromagnetic regime, where < T^j^j^y and the 
heavy fermion regime, where 7^ > Tj^j^j^y- Experiment has told us in recent times that the transition 
between these two regimes is a quantum critical point. The effective Fermi temperature of the heavy 
Fermi liquid is indicated as a solid line. Circumstantial experimental evidence suggests that this scale 
drops to zero at the antiferromagnetic quantum critical point, but this is still a matter of controversy. 



scattering at each site will form a highly renormalized band, of width ~ 7^. By contrast 
to the single impurity Kondo effect, in the heavy electron phase of the Kondo lattice the 
strong elastic scattering at each site acts in a coherent fashion, and does not give rise 
to a resistance. For this reason, as the heavy electron state forms, the resistance of the 
system drops towards zero. One of the fascinating aspects of the Kondo lattice concerns 
the Luttinger sum rule. This aspect was first discussed in detail by Martin[|3^, who 
pointed out that the Kondo model can be regarded as the result of adiabatically increasing 
the interaction strength U in the Anderson model, whilst preserving the valence of the 
magnetic ion. During this process, one expects sum rules to be preserved. In the impurity, 
the scattering phase shift at the Fermi energy counts the number of localized electrons, 
according to the Friedel sum rule 



This sum rule survives to large U, and reappears as the constraint on the scattering 
phase shift created by the Abrikosov Suhl resonance. In the lattice, the corresponding 
sum rule is the Luttinger sum rule, which states that the Fermi surface volume counts 
the number of electrons, which at small U is just the number of localized (4f, 5f or 3d) 
and conduction electrons. When U becomes large, number of localized electrons is now 
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FIGURE 15. Development of coherence in heavy fermion systems. Resistance in Ce^_^LaxCu^ after 
Onuki and Komatsubara[|3^ 



the number of spins, so that 



2—^ = ne + n 



'{IkY ' ' ' "spins 

This sum rule is thought to hold for the Kondo lattice Hamiltonian, independently of the 
origin of the localized moments. Such a sum rule would work, for example, even if the 
spins in the model were derived from nuclear spins, provided the Kondo temperature 
were large enough to guarantee a paramagnetic state. 

Experimentally, there is a great deal of support for the above picture. It is possible, 
for example, to examine the effect of progressively increasing the concentration of Ce 
in the non-magnetic host LaCUf^.<^^ ) At dilute concentrations, the resistivity rises to a 
maximum at low temperatures. At dense concentrations, the resistivity shows the same 
high temperature behavior, but at low temperatures coherence between the sites leads to 
a dramatic drop in the resistivity. The thermodynamics of the dense and dilute system 
are essentially identical, but the transport properties display the effects of coherence. 

There are also indications that the Fermi surface of heavy electron systems does 
have the volume which counts both spins and conduction electrons, derived from Fermi 



surface studies. [52, 53] 



3.1. Some difficulties to overcome. 

The Doniach scenario for heavy fermion development is purely a comparison of en- 
ergy scales: it does not tell us how the heavy fermion phase evolves from the antiferro- 
magnet. There were two early objections to Doniach's idea: 

• Size of the Kondo temperature Tf^. Simple estimates of the value of Jp required for 
heavy electron behavior give a value 7p ~ 1 . Yet in the Anderson model, 7p ~ 1 
would imply a mixed valent situation, with no local moment formation. 



• Exhaustion paradox. The naive picture of the Kondo model imagines that the local 
moment is screened by conduction electrons within an energy range of the Fermi 
energy. The number of conduction electrons in this range is of order Tj^/D << 1 
per unit cell, where D is the band-width of the conduction electrons, suggesting that 
there are not enough conduction electrons to screen the local moments. 

The resolution of these two issues are quite intriguing. 



3.1.1. Enhancement of the Kondo temperature by spin degeneracy 

The resolution of the first issue has its origins in the large spin-orbit coupling of 
the rare earth or actinide ions in heavy electron systems. This protects the orbital 
angular momentum against quenching by the crystal fields. Rare earth and actinide ions 
consequently display a large total angular momentum degeneracy = 2j + 1, which has 
the effect of dramatically enhancing the Kondo temperature. Take for example the case 
of the Cerium ion, where the 4/^ electron is spin-orbit coupled into a state with j = 5/2, 
giving a spin degeneracy ofA'^ = 27-|-l = 6. Ytterbium heavy fermion materials involve 
the Yb : 4/^^ configuration, which has an angular momentum j = 7/2, or = 8. 

To take account of these large spin degeneracies, we need to generalize the Kondo 
model. This was done in the mid-sixties by Coqblin and Schrieffer[p3|]. Coqblin and 
Schrieffer considered a degenerate version of the infinite U Anderson model in which 
the spin component of the electrons runs from — j to j, 

k(j CT k,a 

Here the conduction electron states are also labeled by spin indices that run from — j 
to j. This is because the spin-orbit coupled / states couple to partial wave states of the 
conduction electrons in which the orbital and spin angular momentum are combined into 
a state of definite j. Suppose \ko) represents a plane wave of momentum k, then one can 
construct a state of definite orbital angular momentum / by integrating the plane wave 
with a spherical harmonic, as follows: 

da. 



\klma) = j — \ka)YUk) 



When spin orbit interactions are strong, one must work with a partial wave of definite 
j, obtained by combining these states in the following linear combinations. Thus for the 
case J = Z + 1/2 (relevant for Ytterbium ions), we have 



M= L \i 21+1 ' '^^"^-2-2^- 



tT=±l 



An electron creation operator is constructed in a similar way. This construction is 
unfortunately, not simultaneously possible at more than one site. 



When Ej «Q, the valence of the ion approaches unity and 1. In this limit, 

one can integrate out the virtual fluctuations ^ + e~ via a Schrieffer Wolff 
transformation. This leads to the Coqblin Schrieffer model 

^CS = llh^\o^ko+'f E '^\fh!a^a^^ ((7, a, ^ G [-7,7]). 

ka k,k',afi 

where J = V'^ l\Ej \ is the induced antiferromagnetic interaction strength. This interaction 

is understood as the result of virtual charge fluctuations into the state, ^ + 
The spin indices run from — j to j, and we have introduced the notation 



Notice that the charge Q = n^oi the /—electron, normally taken to be unity, is conserved 
by the spin-exchange interaction in this Hamiltonian. 

To get an idea of how the Kondo effect is modified by the larger degeneracy, consider 
the renormalization of the interaction, which is given by the diagram 




D 



= J + Nrpln\^-j (86) 

( where the cross on the intermediate conduction electron state indicates that all states 
with energy le^^l e [D'^D] are integrate over). From this result, we see that j8(g) = 
dg{D)/d\nD = —Ng^, where g — Jp has an N— fold enhancement, derived from the 
intermediate hole states. A more extensive calculation shows that the beta function to 
third order takes the form 

I5{g) = -Ng^+Ng\ (87) 
This then leads to the Kondo temperature 



Tj^ = D{NJp)^exp 



1 



NJp 



so that large degeneracy enhances the Kondo temperature in the exponential factor. 
By contrast, the RKKY interaction strength is given by Tj^j^^^y ~ fip, and it does not 
involve any A'^ fold enhancement factors, thus in systems with large spin degeneracy, 
the enhancement of the Kondo temperature favors the formation of the heavy fermion 
ground-state. 

In practice, rare-earth ions are exposed to the crystal fields of their host, which splits 
ihc N — 2j -\- I fold degeneracy into many multiplets. Even in this case, the large 
degeneracy is helpful, because the crystal field splitting is small compared with the band- 
width. At energies D' large compared with the crystal field splitting T^, D' » T^, the 
physics is that of an A' fold degenerate ion, whereas at energies D' small compared with 
the crystal field splitting, the physics is typically that of a Kramers doublet, i.e. 



from which we see that at low energy scales, the leading order renormalization of g is 
given by 

> 1 .mn[£,y2lJ^' 



where the first logarithm describes the high energy screening with spin degeneracy A'^, 
and the second logarithm describes the low-energy screening, with spin degeneracy 2. 
This expression is ~ when Z)' ~ T^, the Kondo temperature, so that 



T I \ T* 



from which we deduce that the renormalized Kondo temperature has the form[^ 

1 



= Dexp 



N 

1 \ /Z)\ 2 



2JopJ 



Here the first term is the expression for the Kondo temperature of a spin 1/2 Kondo 
model. The second term captures the enhancement of the Kondo temperature coming 
from the renormalization effects at scales larger than the crystal field splitting. Suppose 
~ 100^, and D ~ 1000^, and N = 6, then the enhancement factor is order 100. This 
effect enhances the Kondo temperature of rare earth heavy fermion systems to values 
that are indeed, up to a hundred times bigger than those in transition metal systems. This 
is the simple reason why heavy fermion behavior is rare in transition metal systems. 



[ pSQ In short- spin-orbit coupling, even in the presence of crystal fields, substantially 



enhances the Kondo temperature. 



3.1.2. The exhaustion problem 

At temperatures T^T^^, a local moment is "screened" by conduction electrons. What 
does this actually mean? The conventional view of the Kondo effect interprets it in terms 
of the formation of a "magnetic screening cloud" around the local moment. According to 
the screening cloud picture, the electrons which magnetically screen each local moment 
are confined within an energy range of order 5e ~ around the Fermi surface, giving 
rise to a spatially extended screening cloud of dimension I = Vp/Tf. ~ a^, where a is a 
lattice constant and £^ is the Fermi temperature. In a typical heavy fermion system, this 



length would extend over hundreds of lattice constants. This leads to the following two 
dilemmas 



1. It suggests that when the density of magnetic ions is greater than p ~ l/P, the 
screening clouds will interfere. Experimentally no such interference is observed, 
and features of single ion Kondo behavior are seen at much higher densities. 

2. " The exhaustion paradox" The number of "screening"electrons per unit cell within 
energy of the Fermi surface roughly Tj^/W , where W is the bandwidth, so 
there would never be enough low energy electrons to screen a dense array of local 
moments. 

In this lecture I shall argue that the screening cloud picture of the Kondo effect 
is conceptually incorrect. Although the Kondo effect does involve a binding of lo- 
cal moments to electrons, the binding process takes place between the local moment 
and high energy electrons , spanning decades of energy from the Kondo temperature 
up to the band-width. (Fig. [1^) I shall argue that the key physics of the Kondo ef- 
fect, both in the dilute impurity and dense Kondo lattice, involves the formation of a 
composite heavy fermion formed by binding electrons on logarithmically large energy 
scales out to the band-width. These new electronic states are injected into the conduction 
electron sea near the Fermi energy. For a single impurity, this leads to a single isolated 
resonance. In the lattice, the presence of a new multiplet of fermionic states at each 
site leads to the formation of a coherent heavy electron band with an expanded Fermi 
surface. ( ITH) 



3.2. Large N Approach 

We shall now solve the Kondo model, both the single impurity and the lattice, in 
the large limit. In the early eighties, Anderson pointed out that the large spin 
degeneracy N = 2j+l furnishes a small parameter 1/A^ which might be used to develop 
a controlled expansion about the limit N ^ o<^. Anderson's observation immediately 
provided a new tool for examining the heavy fermion problem: the so called "large A'^ 
expansion". pPj]. 

The basic idea behind the large expansion, is to take a limit where every term in the 
Hamiltonian grows extensively with A^^. In this limit, quantum fluctuations in intensive 
variables, such as the electron density, become smaller and smaller, scaling as 1/A^, and 
in this sense, 

N 

behaves as an effective Planck's constant for the theory. In this sense, a large A'^ expansion 
is a semi-classical treatment of the quantum mechanics, but instead of expanding around 
H = 0, one can obtain new, non trivial results by expanding around the non trivial solvable 
limit ^ = 0. For the Kondo model, we are lucky, because the important physics of the 
Kondo effect is already captured by the large A^ limit as we shall now see. 
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FIGURE 16. Contrasting (a) the "screening cloud" picture of the Kondo effect with (b) the composite 
fermion picture. In (a), low energy electrons form the Kondo singlet, leading to the exhaustion problem. 
In (b) the composite heavy electron is a highly localized bound-state between local moments and high 
energy electrons which injects new electronic states into the conduction sea at the chemical potential. 
Hybridization of these states with conduction electrons produces a singlet ground-state, forming a Kondo 
resonance in the single impurity, and a coherent heavy electron band in the Kondo lattice. 



Our model for a Kondo lattice or an ensemble of Kondo impurities localized at sites 
7 is 

" = L^/u^U + I."iU) (89) 
ka J 

where 

is the interaction Hamiltonian between the local moment and conduction sea. Here, the 
spin of the local moment at site j is represented using pseudo-fermions 

and 

k 

creates an electron localized at site j. 



There are a number of technical points about this model that need to be discussed: 



• The spherical cow approximation. For simplicity, we assume that electrons have 
a spin degeneracy = 2j + 1. This is a theorists' idealization- a "spherical cow 
approximation" which can only be strictly justified for a single impurity. Neverthe- 
less, the basic properties of this toy model allow us to understand how the Kondo 
effect works in a Kondo lattice. With an A^-fold conduction electron degeneracy, it 
is clear that the Kinetic energy will grow as 0{N). 

• Scaling the interaction. Now the interaction part of the Hamiltonian Hj{j) in- 
volves two sums over the spin variables, giving rise to a contribution that scales as 
0{N^). To ensure that the interaction energy grows extensively with A^, we need to 
scale the coupling constant as 0{\/N). 

• Constraint = Q. Irreducible representations of the rotation group SU (N) re- 
quire that the number of /—electrons at a given site is constrained to equal to 

= Q. In the large limit, it is sufficient to apply this constraint on the aver- 
age (hj:) = Q, though at finite A'^ a time dependent Lagrange multiplier coupled 
to the difference rij: — Qh required to enforce the constraint dynamically. With Q 

/—electrons, the spin operators Y^^^ = f^f^^ provide an irreducible antisymmet- 
ric representation of SU (N) that is described by column Young Tableau with Q 
boxes. As A^ is made large, we need to ensure that q = Q/N remains fixed, so that 
Q ~ 0{N) is an extensive variable. Thus, for instance, if we are interested in A^ = 2, 
this corresponds to q = n^/N = j. We may obtain insight into this case by consid- 
ering the large A^ limit with q = 1/2. 

The next step in the large A'^ limit is to carry out a "Hubbard Stratonovich" transfor- 
mation on the interaction. We first write 



N 

with a summation convention on the spin indices. We now factorize thisp5|, B5p as 



This is an exact transformation, provided the hybridization variables Vj{t) are regarded 
as fluctuating variables inside a path integral, so formally. 



zixy] 



^[y,A]Tr[rexp 







(90) 



where 



//[y, A] = iy^,,c.,, + 1 J] + A^.[.,(j) - Q]) , (91) 

ka J 



is exact. In this expression, ^[V, A] denotes a path integral over all possible time- 
dependences of Vj and A.y(T), and T denotes time ordering. The important point for 
our discussion here however, is that in the large limit, the Hamiltonian entering into 
this path integral grows extensively with A^, so that we may write the partition function 
in the form 



^[y,A]Tr[rexp 



.1^ 

-N / Jf[V,X] 



(92) 



where J^[V,X] — jjH[V,X] ~ 0(1) is an intensive variable in N. The appearance of a 
large factor in the exponential means that this path integral becomes dominated by its 
saddle points in the large limit- i.e, if we choose 

where the saddle point values Vg and Xo are chosen so that 



ainZ[y,A] 



dV 

then in the large A^ limit. 



dlnZ[V,X] 



Vj=Vo,Xj=Xo 



= 



2 = ^^^-^H^oM 



In this way, we have converted the problem to a mean-field theory, in which the fluctuat- 
ing variables Vy(T) and X^{x) are replaced by their saddle-point values. Our mean-field 
Hamiltonian is then 



^r^ k ka ka 

ka 



where n is the number of sites in the lattice. We shall now illustrate the use of this mean- 
field theory in two cases- the Kondo impurity, and the Kondo lattice. In the former, there 
is just one site; in the latter, translational invariance permits us to set Vj — Vo at every 
site, and for convenience we shall choose this value to be real. 



3.3. Mean-field theory of the Kondo impurity 

3.3.1. Diagonalization ofMF Hamiltonian 

The Kondo effect is at heart, the formation of a many body resonance. To understand 
this phenomenon at its conceptually simplest, we begin with the impurity model. We 
shall begin by writing down the mean-field Hamiltonian for a single Kondo ion 

H = L^k^\.^ka + Lnc\ja + f\c,^]+^'Lnf^-^Q+ ^ (93) 
ka ka a 



By making a mean-field approximation, we have reduced tlie problem to one of a 
self-consistently determined resonant level model. Now, suppose we diagonalize this 
Hamiltonian, writing it in the form 



H = ^ EyO^ yfjUyfj H " XQ 



(94) 



where the "quasiparticle operators" Uj are related via a unitary transformation to the 
original operators 

a\a = Y.^/k<,+^f^o- (95) 

k 



commuting y^j with H, we obtain 



(96) 



Expanding the right and left-hand side of ( |9q ) in terms of ( |95| ) and (|93|), we obtain. 







-y£a, + (£^-A)/3 = 



(97) 



Solving for a^. using the first equation, and substituting into the second equation, we 
obtain 

^7-^-1^^ = (98) 

k H 

We could have equally well obtained these eigenvalue equations by noting the electron 



eigenvalues Ej must correspond to the poles of the f-Green function, Gj:{Ey) 
where from an earlier subsection. 



0, 



CO 



k 



(99) 



Either way, the one-particle excitation energies Ey must satisfy 



y2 



k ^k 



(100) 



The solutions of this eigenvalue equation are illustrated graphically in Fig. ([T7|). Suppose 
the energies of the conduction sea are given by the 2M discrete values 

£j^ = {k+l-)A£, kE{-M,...,M-l} 



Suppose we restrict our attention to the particle-hole case when the f-state is exactly 
half filled, i.e. when Q = N/2. In this situation, A = 0. We see that one solution to the 
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FIGURE 17. (a)Graphical solution of the equation y = — -^t^ for eight equally spaced conduction 
electron energies. Notice how the introduction of a new bound-state at y = displaces electron band-states 
both up and down in energy. In this way, the Kondo effect injects new bound-state fermion states into the 
conduction sea. (b) Energy dependence of the scattering phase shift. 

eigenvalue equation corresponds to Ey = 0. The original band-electron energies are now 
displaced to both lower and higher energies, forming a band of 2M + 1 eigenvalues. 
Clearly, the effect of the hybridization is to inject one new fermionic eigenstate into the 
band. Notice however, that the electron states are displaced symmetrically either-side of 
the new bound-state at Ey = 0. 



Each new eigenvalue is shifted relative to the original conduction electron energy by 
an amount of order Ae. Let us write 



5y 

Ey = £y — Ac 

' ' K 

where 5 e [0, it] is called the "phase shift". Substituting this into the eigenvalue equation, 
we obtain 

Ey = A + ^ 

n=y+l-M — ^) 

Now if M is large, we can replace the sum over states in the above equation by an 
unbounded sum 

V} ^ 1 



' A8„i^^(„-|) 
Using contour integration methods, one can readily show that 



— ;rcot5 



so that the phase shift is given by 5y = d{Ey), where 

KpV^ 



tan5[e] 



where we have replaced p = ^ as the density of conduction electron states. This can 
also be written 

A 



5(e) = tan ^ 



A -e 



= Imln[A + jA-e] (101) 



where A = npV^ is the width of the resonant level induced by the Kondo effect. Notice 
that for A = 0, 5 = ;r/2 at the Fermi energy. 

• The phase shift varies from 5 = at = — oo to 5 = ;r at = oo , passing through 
5 = ;r/2 at the Fermi energy. 

• An extra state has been inserted into the band, squeezing the original electron states 
both down and up in energy to accommodate the additional state: states beneath 
the Fermi sea are pushed downwards, whereas states above the Fermi energy are 
pushed upwards. From the relation 

Ey S-y 5 {^Ey^ 



we deduce that 
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= ^ + ——7^ (102) 
np dE 

where p = 1 / Ae is the density of states in the continuum. The new density of states 
p*(£')is given by p*{E)dE = pde, so that 

P*iE)=piO)^ = P + P,iE) (103) 

where 

, , 1 dd(E) 1 A 

corresponds to the enhancement of the conduction electron density of states due to 
injection of resonant bound-state. 



3.3.2. Minimization of Free energy 

With these results, let us now calculate the Free energy and minimize it to self- 
consistently evaluate A and A. The Free energy is given by 

F = -NTY,H^+e-^^y]-XQ + ^. (105) 
In the continuum limit, where e ^ 0, we can use the relation E'y = Cy — Aef to write 



^^0 



/ ^ Ae 

= -T\n[l+e-^''y] 5(£y)/(ey) (106) 

K 

where f{x) = 1 / {e^^ + 1 ) is the Fermi function. The first term in ( |106| ) is the Free energy 
associated with a state in the continuum. The second term results from the displacement 
of continuum states due to the injection of a resonance into the continuum. Inserting this 
result into ( |105D , we obtain 

Ae NV'^ 

- Fo-N ^/(e)5(e)-Ae+^ 

The shift in the Free energy due to the Kondo effect is then 

f'°° de NA 
AF = -N — /(e)Imln[C-e]-Ae+— — (108) 
j-oo 71 nJp 
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where we have introduced ^ = A + iA. This integral can be done at finite temperature, 
but for simplicity, let us carry it out at T = 0, when the Fermi function is just at step 
function, f{x) = 6{—x). This gives 
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where we have expanded + D) In 
We can further simplify this expression by noting that 



Din [ + ^ InD to obtain the second line. 
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where q = Q/N. With this simplification, the shift in the ground-state energy due to the 
Kondo effect is 



AE = -Im 



Cln 



(111) 



where we have dropped the constant term and introduced the Kondo temperature 
De^ . The stationary point dE/dC, = is given by 



C = A+/A 



ta.n{nq) 



VA2 + A2 
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Notice that 



The phase shift 5 = nq is the same in each spin scattering channel, reflecting 
the singlet nature of the ground state. The relationship between the filling of the 
resonance and the phase shift Q = Y^^y = is nothing more than Friedel's sum 
rule. 

The energy is stationary with respect to small variations in A and A. It is only a local 
minimum once the condition dE / dX, corresponding to the constraint (hj:) = Q, or 
A = Acot{7i:q) is imposed. It is instructive to study the energy for the special case 
q — ^, X — which is physically closest to the S = 1/2, N = 2 case. In this case, 
the energy takes the simplified form 



A In 



eT, 



(112) 



Plotted as a function of V , this is the classic "Mexican Hat" potential, with a 
minimum where dE / dV = OatA = ;rp|yp = T^. (Fig. [IS]) 

According to ( |103D , the enhancement of the density of states at the Fermi energy is 



= p + 
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per spin channel. When the temperature is changed or a magnetic field introduced, 
one can neglect changes in A and A, since the Free energy is stationary. This implies 
that in the large N limit, the susceptibility and linear specific heat are those of a non- 
interacting resonance of width A. The change in linear specific heat ACy = Ayr and 
the change in the paramagnetic susceptibility Ax are given by 
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Notice how it is the Kondo temperature that determines the size of these two 
quantities. The dimensionless "Wilson" ratio of these two quantities is 



W = 



Ay 



At finite A'^, fluctuations in the mean-field theory can no longer be ignored. These 
fluctuations induce interactions amongst the quasiparticles, and the Wilson ratio 
becomes 

1 

W 



1 



I' 

N 



The dimensionless Wilson ratio of a large variety of heavy electron materials lies 
remarkably close to this value. 



3.4. Gauge invariance and the composite nature of the /—electron 

We now discuss the nature of the /—electron. In particular, we shall discuss how 

• the /—electron is actually a composite object, formed from the binding of high- 
energy conduction electrons to the local moment. 

• although the broken symmetry associated with the large mean-field theory does 

not persist to finite N, the phase stiffness associated with the mean-field theory 
continues to finite A'^. This phase stiffness is responsible for the charge of the 
composite / electron. 



3.4.1. Composite nature of the heavy f— electron 

Let us begin by discussing the composite structure of the /—electron. In real mate- 
rials, the Kondo effect we have described involves spins formed from localized f- or 
d-electrons. Though it is tempting to associate the composite /—electron in the Kondo 



effect with the the /—electron locked inside the local moment, we should also bear in 
mind that the Kondo effect could have occurred equally well with a nuclear spin! Nu- 
clear spins do couple antiferromagnetically with a conduction electron, but the coupling 
is far too small for an observable nuclear Kondo effect. Nevertheless, we could conduct a 
thought experiment where a nuclear spin is coupled to conduction electrons via a strong 
antiferromagnetic coupling. In this case, a resonant bound-state would also form from 
the nuclear spin. The composite bound-state formed in the Kondo effect clearly does not 
depend on the origin of the spin partaking in the Kondo effect. 

There are some useful analogies between the formation of the composite /—electron 
in the Kondo problem and the formation of Cooper pairs in superconductivity, which 
we shall try to draw upon. One of the best examples of a composite bound-state is the 
Cooper pair. Inside a superconductor, pairs of electrons behave as composite bosonic 
particles. One of the signatures of pair formation, is the fact that Cooper pairs of electron 
operators behave as a single composite at low energies, 

VAj,(x)v^^(x') = F{x-x') 

The Cooper pair operator is a boson, and it behaves as a c-number because the Cooper 
pairs condense. The Cooper pair wavefunction is extremely extended in space, extending 
out to distances of order ^ ~ Vp/Tc. A similar phenomenon takes place in the Kondo 
effect, but here the bound-state is a fermion and it does not condense For the Kondo 
effect the fermionic composite {a • S(x))^^ va^(x) behaves as a single charged electron 
operator. The analogy between superconductivity and the Kondo effect involves the 
temporal correlation between spin-flips of the conduction sea and spin-flips of the local 
moment, so that at low energies 

[o„p-Sit)]xifpit')r. Ait -t') fait'). 

The function A{t — t') is the analog of the Cooper pair wavefunction, and it extends out 
to times r^H/T^. 

To see this in a more detailed fashion, consider how the interaction term behaves. In 
the path integral we factorize the interaction as follows 

By comparing these two terms, we see that the composite operator ^f^p{j)YaU) be- 
haves as a single fermi field: 



Evidently, a localized conduction electron is bound to a spin-Hip of the local moment 
at the same site, creating a new independent fermionic excitation. The correlated action 
of adding a conduction electron with a simultaneous spin flip of the local moment at 
the same site creates a composite f— electron. 



It is worth noting that this fermionic object only hybridizes with conduction electrons at 
a single point: it is thus local in space. 

Let us now try to decompose the composite fermion in terms of the electrons that 
contribute to the bound-state amplitude. We start by writing the local moment in the 
fermionic representation, ^ 



where we have replaced the bilinear product between the conduction and /—electron 
by its expectation value. We can evaluate this "bound-state amplitude" from the corre- 
sponding Green-function 
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where we have chosen the half-filled case Q/N = 1/2, A = 0. In the large band-width 
limit, the main contribution to this integral is obtained by neglecting the principal part 
of the conduction electron propagator l/(ft) — e, — i5) — > in5{co — e, ), so that 



(116) 



From this expression, we can see that the contribution of a given k state in the Fermi sea 
to the bound-state amplitude is given by 



-(/' 



N 



e2 + A2 



This function decays with the inverse of the energy, right out to the band- width. Indeed, 
if we break-down the contribution to the overall bound-state amplitude, we see that each 
decade of energy counts equally. Let us take 7 = and divide the band on a logarithmic 
scale into n equal parts, where the ratio of the lower and upper energies is 5 > 1, then 
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^ Important and subtle point: The emergence of a composite fermion does not depend on a fermionic 
representation of the spin. The fermionic representation for the spin is simply the most convenient because 
it naturally furnishes us with an operator in the theory that represents the composite bound-state. 



r Ds " 1 

= pVq^^ ln5 + lns+...lns + ln^— > (117) 

This demonstrates that the composite bound-state involves electrons spread out over 
decades of energy out to the band-width. If we complete the integral, we find that 

as expected from the minimization of the energy. Another way of presenting this discus- 
sion, is to write the composite bound-state in the time-domain, as 



where now 



],T^p{t')xifa(t)^Kt-t')fa{t') (118) 



^{t-t') = ^{f;{t)w^{t')) 



This is the direct analog of Cooper pair bound-state wavefunction, except that the 
relevant variable is time, rather than space. If one evaluates the function A(?) at a finite 
t, we find that 



k \^k 



Heuristically, the finite time cuts off the energy integral over the Fermi surface at an 
energy of ordexH/t, so that 



/ P^oln(f) {t«n/T^) 
^(^)'^\pV.ln(^) {t»n/T^) 



emphasizing the fact that the Kondo effect involves a correlation between the spin-flips 
of the conduction sea and the local moment over decades of time scales from the the 
inverse band- width up to the Kondo time^/r^. 

From these discussions, we see that the Kondo effect is 

• entirely localized in space. 

• extremely non-local in time and energy. 

This picture of the Kondo effect as a temporal, rather than a spatial bound-state is vital 
if we are to understand the extension of the Kondo effect from the single impurity to the 
lattice. 



3.4.2. Gauge invariance and the charge of the f— electron 

One of the interesting points to emerge from the mean-field theory is that the energy of 
mean-field theory does not depend on the phase of the bound-state amplitude V = |y 



This is analogous to the gauge invariance in superconductivity, which derives from 
the conservation of the total electronic charge. Here, gauge invariance arises because 
there are no charge fluctuations at the site of the local moment, a fact encoded by the 
conservation of the total f-charge Q. Let us look at the full Lagrangian for the /—electron 
and interaction term 

= faHidt-^)fa-Hj 
Hi = v(^xt/\fa) + [faVa)v+N^ (119) 



This is invariant under the "Read-Newns"[|4g] transformation 



V ^ ye'^ {e^G + ^), 

X A + (120) 
ot 

where the last relation arises from a consideration of the gauge invariance of the dynamic 
part f\idj - A)/ of the Lagrangian. Now if V{t) = \V{t)\e'^^'\ where V{t) is real. 
Read and Newns observed that by making the gauge choice ^{t) = —6{t), the resulting 
V = |y|e'(e+^) = \V\ is real. In this way, once the Kondo effect takes place the phase 
of V = |y is dynamically absorbed into the constraint field A : effectively A = (9^0 
represents the phase precession rate of the hybridization field. The absorption of the 
phase of an order parameter into a dynamical gauge field is called the "Anderson 
Higgs " mechanism. IgTp By this mechanism, once the Kondo effect takes place, V 
behaves as a real, and hence neutral object under gauge transformations, this in turn 
implies that the composite /—electron has to transform under real electromagnetic gauge 
transformations, in other words the Anderson Higgs effect in the Kondo problem endows 
the composite /—electron with charge. 

There is a paradox here, for in the Kondo effect, there can actually be no true broken 
symmetry, since we are dealing with a system where the number of local degrees of free- 
dom is finite. Nevertheless , the phase does develop a stiffness- a stiffness against vari- 
ation in time, and the order parameter consequently develops infinite range correlations 
in time. There is a direct analogy between the spatial phase stiffness of a superconductor 
and the temporal phase stiffness in the Kondo effect. In superconductivity, the energy 
depends on spatial derivatives of the phase 

D 1 2 

£-^(V<^)-2M)2^— ocp, 
2 A^ 

( where we have set^ =1.) Gauge invariance links this stiffness to the mass of the 
photon field, which generates the Meissner effect; the inverse squared penetration depth 
is directly proportional to the phase stiffness. In an analogous fashion, in the Kondo 
effect, the energy depends on temporal derivatives of the phase and the phase stiffness 




FIGURE 18. "Mexican Hat Potential" which determines minimum of Free energy, and self-consistently 
determines the width of the Kondo resonance. The Free energy displays this form provided the constraint 
dF/dX = («^) — 2 = is imposed. 



IS^ 

For a Kondo lattice, there is one independent Kondo phase for each spin site, and 
the independent conservation of Q at each site guarantees that there is no spatial phase 
stiffness associated with 0. The temporal phase stiffness leads to a slow logarithmic 
growth in the phase -phase correlation functions, which in turn leads to power-law 
temporal correlations in the order parameter V(t): 

(5(|)(t)50(t')) ~^Mt-t'), (y(T)y(TO)~e"<^'^(^^^'^(^')^ ~(T-T')~i 

In this respect, the Kondo ground-state resembles a two dimensional superconductor, 
or a one dimensional metal: it is critical but has no true long-range order. As in the 



Note that because A ~ d,^, the phase stiffness is given by = d^F/dX^ 



TABLE 1. Parallels between Superconductivity and the Kondo effect . 
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superconductor, the development of phase stiffness involves real physics. When we make 
a gauge transformation of the electromagnetic field, 

e<i>{x,t) e^{x,t) +dta{x,t), 
eA{x,t) eA{x,t) + Va, 



(121) 



Because of the Anderson - Higgs effect, the hybridization is real and the only way 
to keep Lj invariant under the above transformation, is by gauge transforming the 
/—electron and the constraint field 

A ^ X+dta (122) 
( Notice how A transforms in exactly the same way as the potential e4>.) 

The non-trivial transformation of the /—electron under electromagnetic gauge trans- 
formations confirm that it has acquired a charge. Rigidity of the Kondo phase is thus 
intimately related to the formation of a composite charged fermion. The gauge invariant 
form for the energy dependence of the Kondo effect on the Kondo phase ^ must then be 

£oc^(a,0-e4>)2 

From the coefficient of 4>^, we see that the Kondo cloud has an intrinsic capacitance 
C — e^p^ {E ~ C4>^/2). But since the energy can also be written (eny)^/2C ~ U*nj-/2 
we see that the stiffness of the Kondo phase can also be associated with an interaction 
between the /—electrons of strength U*, where 

3.5. Mean-field theory of the Kondo Lattice 

3.5.1. Diagonalization of the Hamiltonian 

We can now make the bold jump from the single impurity problem, to the lattice. 
Most of the methods described in the last subsection generalize very naturally from the 
impurity to the lattice: the main difficulty is to understand the underlying physics. The 
mean-field Hamiltonian for the lattice[^, ^] takes the form 

where is the number of sites in the lattice. Notice, before we begin, that the composite 
f-state at each site of the lattice is entirely local, in that hybridization occurs at one site 
only. Were the composite f-state to be in any way non-local, we would expect that the 
hybridization of one f-state would involve conduction electrons at different sites. We 
begin by rewriting the mean field Hamiltonian in momentum space, as follows 




where 
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is the Fourier transform of the /—electron field. The absence of k— dependence in 
the hybridization is evident that each composite /—electron is spatially local. This 
Hamiltonian can be diagonalized in the form 
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where a^- and b^^ are linear combinations of ct and , playing the role of 

ka ka ka ■' ka^ ^ ^ 

"quasiparticle operators" of the theory and the momentum state eigenvalues E^_^ of this 
Hamiltonian are determined by the condition 
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which gives 
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are the energies of the upper and lower bands. The dispersion described by these energies 
is shown in Fig. . A number of points can be made about this dispersion: 



We see that the Kondo effect injects new fermionic states into the original con- 
duction band. Hybridization between the heavy electron states and the conduction 
electrons builds an upper and lower Fermi band separated by a "hybridization gap" 
of width Ag = Eg{+) — Eg{—), such that energies in the range 



Eg{-)< E <Xo+Eg{+) 
Egi±) = Xo±^ 



(124) 



are forbidden. Here are the top and bottom of the conduction band. In the 
special case where Xq = 0, corresponding to half filling, a Kondo insulator is 
formed. 

The effective mass of the Fermi surface has the opposite sign to the original con- 
duction sea from which it is built, so naively, the Hall constant should change sign 
when coherence develops. 

The Fermi surface volume expands in response to the presence of the new heavy 
electron bands. The new Fermi surface volume now counts the total number of 
particles. To see this note that 



kla 




FIGURE 19. (a) Dispersion produced by the injection of a composite fermion into the conduction sea. 
(b) RenormaUzed density of states, showing "hybridization gap" (Ag). (c) Transformation of the Fermi 
surface from a light electron Fermi surface into a heavy "hole"-like Fermi surface. 



where n. 



''kXc ~ ^ kXc^kXc number operator for the quasiparticles and ric is the 
total number of conduction electrons. This means 



Ntot = N 
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= Q + nc 



This expansion of the Fermi surface is a direct manifestation of the creation of 
new states by the Kondo effect. It is perhaps worth stressing that these new states 
would form, even if the local moments were nuclear in origin. In other words, it is 
only the rotational degrees of freedom of the local moments that are needed to form 
heavy electron bound-states with the conduction electrons. 

The Free energy of this system is then 
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Let us discuss the ground-state energy. Eg -the T —>^0 Umit of this expression. We can 
write this in the form 



dEp*{E)E+(^-Xq^ 



where we have introduced the density of heavy electron states p*{E) = 52^_^5(£' — 
e[^^). Now the relationship between the energy of the heavy electrons (E) and the 

k 

energy of the conduction electrons (e) is given by 

E = e + 



E-X 



so that the density of heavy electron states related to the conduction electron density of 
states p by 

+ (125) 

The originally flat conduction electron density of states is now replaced by a "hy- 
bridization gap", flanked by two sharp peaks of width approximately KpV^ ~ 7^. With 
this information, we can carry out the integral over the energies, to obtain 
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where we have assumed that the upper band is empty, and the lower band is partially 

_ j_ 

filled, and set Tj^ = De Jp as before. If we impose the constraint 

If = 2 = (127) 

we obtain 

A 

^-''=» 

so that the ground-state energy can be written 

En A / Ae , 

-^ = -ln — ^ . (128) 



Let us pause for a moment to consider this energy functional qualitatively. The Free 

energy surface has the form of the "Mexican Hat" at low temperatures. The minimum 
of this functional will then determine a family of saddle point values V = Voe^^, where 



6 can have any value. If we differentiate the ground-state energy with respect to V^, we 
obtain 

^_ 1 Ae2 \ 

or 

A = ^Tj, (130) 

confirming that A ~ T^. 



3.5.2. Composite Nature of the heavy quasiparticle in the Kondo lattice. 



We now turn to discuss the nature of the heavy quasiparticles in the Kondo lattice. 
Clearly, at an operational level, the composite /—electrons are formed in the same way 
as in the impurity model, but at each site, i.e 



1 



This composite object admixes with conduction electrons at a single site- site j. The 
bound-state amplitude in this expression can be written 
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(131) 



To evaluate the contributions to this sum, it is useful to notice that the condition 
dE/dV = can be written 
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where we have used ( |126p to evaluate the derivative. From this we see that we can write 
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It is clear that as in the impurity, the composite /—electrons in the Kondo lattice 
are formed from high energy electron states all the way out to the bandwidth. In a 
similar fashion to the impurity, each decade of energy between and D contributes 
equally to the overall bound-state amplitude. The above expression only differs from the 
corresponding impurity expression ( |1 15| ) at low energies, showing that low energy elec- 
trons play a comparatively unimportant role in forming the composite heavy electron. It 



is this feature that permits a dense array of composite fermions to co-exist throughout 
the crystal lattice. 

These composite /—electrons admix with the conduction electrons to produce a heavy 
electron band with a density of states given by ( |125D , 

... d£ / 



dE ^ ' {E-IY 
which, setting E = and using (|127D and (|130D , becomes 

K 

at the Fermi energy. The mass enhancement of the heavy electrons is then 

m* qe^ qD 

— 1 H ~ — 

m pTj^ 

This large factor in the effective mass enhancement can be as much as 1000 in the most 
severely renormalized heavy electron systems. 



3.5.3. Consequences of mass renormalization 

The effective mass enhancement of heavy electrons can be directly observed in a wide 
range of experimental quantities including 



The large renormalization of the linear specific heat coefficient 7* ~ ^7 and Pauli 
susceptibility x* ~ ^X- 

The quadratic temperature (" A" ) coefficient of the resistivity. At low temperatures 
the resistivity of a Fermi liquid has a quadratic temperature dependence, p ~ 

po +AT'^, where A ~ ~ (^) ~ related to the density of three-particle 

excitations. The approximate constancy of the ratio A/ 7^ in heavy fermion systems 



is known as the "Kadowaki-Woods" relation. [pOD 

• The renormalization of the effective mass as measured by dHvA measurements of 
heavy electron Fermi surfaces. [|5T|, |52| , |53| ] 

• The appearance of a heavy quasiparticle Drude feature in the frequency dependent 
optical conductivity a{(o). (See discussion below). 

The optical conductivity of heavy fermion metals deserves special discussion. Ac- 
cording to the f-sum rule, the total integrated optical conductivity is determined by the 
plasma frequency 

dco . , - n Z' ne^ 




where n is the density of electrons. [^In the absence of local moments, this is the total 
spectral weight inside the Drude peak of the optical conductivity. 

What happens to this spectral weight when the heavy electron fluid forms? Whilst we 
expect this sum rule to be preserved, we also expect a new "quasiparticle" Drude peak 
to form in which 

/ dcoa{(o) = f2-— = fi — 
J 2 m* m* 

In other words, we expect the total spectral weight to divide up into a tiny "heavy 
fermion" Drude peak, of total weight /j, where 

ne^ 1 
(y((0) = — --| 



The f-sum rule is a statement about the instantaneous, or short-time diamagnetic response of the metal. 
At short times d j/dt = [ne^ /m)E, so the high frequency limit of the conductivity is o{(0) = ^ But 
using the Kramers Kronig relation 

, , f dx aix) 

a{(o) = / 

J in X — CO — 10 

at large frequencies, 

a{co) = - — — / —(y{x) 
5-1(0 J n 

so that the short-time diamagnetic response implies the f-sum rule. 



is split off by an energy of order V ~ \fT^ from an "inter-band" component associated 
with excitations between the lower and upper Kondo bands. |FR This second term 



carries the bulk ~ of the spectral weight. (Fig. pO|). 

Simple calculations, based on the Kubo formula confirm this basic expectation, 
p4| ] showing that the relationship between the original relaxation rate of the conduction 
sea and the heavy electron relaxation rate T* is 

(T*)-i = £^(T)-i. (134) 

Notice that this means that the residual resistivity 

m* m 

is unaffected by the effects of mass renormalization. This can be understood by observ- 
ing that the heavy electron Fermi velocity is also renormalized by the effective mass, 
= so that the mean-free path of the heavy electron quasiparticles is unaffected by 
the Kondo effect. 

This is yet one more reminder that the Kondo effect is local in space, yet non-local in 
time. 

These basic features- the formation of a narrow Drude peak, and the presence 
of a hybridization gap, have been seen in optical measurements on heavy electron 
systems[|5[ ^ ^ 



3.6. Summary 

In this lecture we have presented Doniach's argument that the enhancement of the 
Kondo temperature over and above the characteristic RKKY magnetic interaction energy 
between spins leads to the formation of a heavy electron ground-state. This enhance- 
ment is thought to be generated by the large spin degeneracies of rare earth, or actinide 
ions. A simple mean-field theory of the Kondo model and Kondo lattice, which ignores 
the RKKY interactions, provides a unified picture of heavy electrons and the Kondo 
effect. The essential physics involves composite quasiparticle formation between high 
energy conduction band electrons and local moments. This basic physical effect is local 
in space, but non-local in time. Certain analogies can be struck between Cooper pair for- 
mation, and the formation of the heavy electron bound-state, in particular, the charge on 
the /—electron can be seen as a direct consequence of the temporal phase stiffness of the 
Kondo bound-state. This bound-state hybridizes with conduction electrons- producing a 
single isolated resonance in a Kondo impurity, and an entire renormalized Fermi surface 
in the Kondo lattice. 



3.7. Exercises 



1. (a) Directly confirm the Read-Newns gauge transformation (120) 



(b) Directly calculate the "phase stiffness" = —j^ of the large Kondo model and 
show that at r = 0. 

N ( sin(7r^)' 



2. (a) Introduce a simple relaxation time into the conduction electron propagator, writing 

i(0„ — X 



G{k,i(On) =i(0„ + hgn{(0„)/2z+- r- (135) 



Show that the poles of this Greens function occur at 

It* 

where 

* m* 

T = T 

m 

is the renormalized elastic scattering time, 
(b) The Kubo formula for the optical conductivity of an isotropic one-band system is 

, , A^^'^ n(v) 

where we have used the fold spin degeneracy, and n(v) is the analytic extension 

n(/v„) = TJ^ G{k, i(0,n) G(k, i(0,n + iVn) - G{k, i(0,n) 



where in our case, G{k,i(On) is the conduction electron propagator. Using (135), and 
approximating the momentum sum by an integral over energy, show that the low 
frequency conductivity of the large Kondo lattice is given by 

cy(v) 



m* (t*)-^ -iv' 

4. QUANTUM CRITICALITY IN HEAVY ELECTRON SYSTEMS 

4.1. Introduction 

This section provides a brief introduction to the unsolved problem of quantum crit- 
ical! ty in heavy fermion materials. Many of the ideas summarized here are the result 
of collaborations, and much of the material in this section is published in review form. 
[ |67i |106| , |107D Heavy electron materials lie on the verge of magnetic instability. In the 



discussion of the last section, we ignored magnetism and focussed on the dense Kondo 
effect. What happens when they are pushed to the very edge of magnetic instability? 
Such a question was first posed in the context of itinerant magnetic order in a pioneering 
work by John Hertz, almost thirty years ago. Hertz concluded that a metallic system at 



the edge of magnetic instability would develop a new kind of critical behavior- quantum 
critical behavior. 

A quantum critical point (QCP) is a zero-temperature instability between two phases 
of matter where quantum fluctuations develop long range correlations in both space and 



time 05 80 ■ At a finite temperature critical point, the critical long-wavelength fluctuations 
of the order parameter do not involve quantum mechanics. This is because thermal 
fluctuations destroy the coherence of quantum fluctuations on time-scales longer than 

T~7^, (136) 

The great revolution in our understanding of critical phenomena which occured in the 
1970s involved many tools borrowed from relativistic field theory, but the physics was 
entirely classical. 

Experimental developments of the past decade have brought a new awareness of the 
importance of quantum critical points in condensed matter physics. These special points 
exert a profound influence on the finite temperature properties of a material. Materi- 
als close to quantum criticality develop a new excitation structure, they display novel 
thermodynamic, transport and magnetic behavior. They also have marked a predeliction 
towards the development of new kinds of order, such as anisotropic superconductivity. 
A dramatic example is provided by the cuprate superconductors. By doping with holes, 
these materials pass through one or more quantum phase transitions: from an insulator 
to a metal with a spin gap at low doping, and at higher doping a second QPT appears 



to occur when the spin gap closes [ p9| ] (Fig. gl]) The singular interactions induced by 
quantum criticality are thought to the driving force for both the high temperature super- 
conductivity and the anomalous metallic state above the spin gap temperature r*.06O|] 

Heavy Fermion materials offer a unique opportunity to study quantum criticality 
under controlled conditions. By the application of pressure, doping and most recently, 
magnetic field, these materials can be tuned through a quantum critical point from a 



metallic antiferromagnet into a paramagnet (Fig. Unlike the cuprate metals, here 
the paramagnetic phase is a well characterized Fermi liquid, 06T|, |62| , [63| ] with heavy 
Landau quasiparticles, or "heavy electrons". A central property of these quasiparticles, 
is the existence of a finite overlap "Z" between a single quasiparticle state, denoted 
by |qp^) and the state formed by adding a single electron to the ground- state, denoted 
by \e^) = c\^\0). This quantity is closely related to the ratio m/m* of the electron to 
quasiparticle mass, 

Z=|(e-|qp-)|2~4. (137) 
m 

A wide body of evidence suggests that m*/m diverges at a heavy fermion QCP, indicat- 
ing that 

Z^O (P^Pc). 
The state which forms above the QCP is referred to as a "non-Fermi" or "singular Fermi 



liquid". [p4 651 57, 58, pq] By what mechanism does this break-down in the Landau 



quasiparticle occur? 



Table. 1. Selected Heavy Fermion compounds with quantum critical points. 




QCP 

FIGURE 21. Schematic phase diagram for cuprate superconductors showing location of possible quan- 
tum critical points. One of these QCP may be responsible for the anomalous normal state which develops 
above the pseudogap scale. 
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FIGURE 22. Illustration of quantum critical physics in heavy fermion metals. As criticality is ap- 
proached from either side of the transition, the temperature scale Tp on which Fermi liquid behavior 
breaks down goes to zero. A key challenge is to characterize the new class of universal excitations which 
develops above Tp . 



New data[78, 75] show a stronger divergence at lower temperatures. 



4.2. Properties of the Heavy Fermion Quantum Critical Point 



There is a growing list of heavy fermion systems that have been tuned to an antifer- 
romagnetic QCP by the application of pressure or by doping (Table 1.). These materials 
display many common properties 

• Fermi liquid behavior in the paramagnet, as indicated by the emergence of a 
quadratic temperature dependence in the resistivity in the approach to the QPT 
p = p^ +AT^ [|8|, ^ at ever lower temperatures. 



Divergent A coefficient in resistivity at the QCP. In a typical Fermi liquid the A 
coefficient in the resistivity is proportional to (^jjf^^ ~ (^) ' ^^^^^ '•^^ 



Fermi temperature. Support for the divergence of the effective mass is provided 
by the observation that the quadratic coefficient A of the resistivity grows, and 
apparantly diverges at the quantum critical point[[79|]. 

Divergent specific heat at the QCP, with an asymptotic logarithmic temperature 
dependence, 



To 



+ %, (138) 



where R is the gas constant, and experimentally Q ^ 0.4, suggesting that the Fermi 
temperature vanishes and the quasiparticle effective masses diverge 

in 

r;^o, ^oo (139) 

m 

at the QCP. The above expression has been written in a form where the character- 
istic energy To enters both inside the logarithm and in the prefactor. There are a 
number of materials where this one-parameter form holds, with 7„ = 0, suggest- 
ing a new kind of universality where no normal component to the Fermi surface 
survives at the QCP. ^] 
Quasi-linear resistivity 

p oc (140) 

at the QCP with £ in the range — 0.6. In critical YbRh2Si2_fiex, p over three 
decades!^ 



Non-Curie spin susceptibilities 

X'\t)=Xo'+cT'' (141) 

with a < 1 observed in critical CeCU(^_^Ux (x=0.1), YbRh2Si2_^Gex (x=0.1) and 
CeNi2Ge2- 

E/T and H/T Scaling. In critical CeCu^ and YbRh2Si2_^Gex the differential 
magnetic susceptibility dM/dH exhibits ///r scaling, 

[dM/dH)-' = Xo' +cT'^g[H/T], (142) 

where a ~ 0.75. Neutron measurements [[7^ show E/T scaling[P9|, in the dy- 
namical spin susceptiblity of critical CeCu^^Aux, throughout the Brillouin zone, 
parameterized in the form 

X-'iq. CO) = T"f{E/T) + Xo'id) (143) 

F[x] oc (1 — ix)". Scaling behavior with a single anomalous exponent in the 
momentum-independent component of the dynamical spin susceptibility suggests 
an emergence of local magnetic moments which are critically correlated in time at 
the quantum critical point [[7^. 



4.3. Universality 



Usually, the physics of a metal above its Fermi temperature depends on the detailed 
chemistry and band- structure of the material: it is non-universal. However, if the renor- 
malized Fermi temperature Tp (P) can be tuned to become arbitrarily small compared 
with the characteristic scales of the material as one approaches a QCP, we expect that 
the "high energy" physics above the Fermi temperature Tp is itself, universal . 

Quantum critical behavior implies a divergence of the long distance and long-time 
correlations in the material. Finite temperatures introduce the cutoff timescale 



H 



T:T = -r^ (144) 



beyond which coherent quantum processes are dephased by thermal fluctuations. Renor- 
malization group principles [PH] imply that the quantum critical physics has an upper- 



critical dimension du- For d < du, Xj becomes the correlation time T of the system[|920, 
so frequency dependent correlation functions and response functions take the form 

F(co,r) = -^noJTp) = l^f{ncD/k^T). (145) 

leading to E/T scaling []93|]. By contrast, for d > du the correlation time is sensitive to 
the details of the short-distance interactions between the critical modes, and in general 
T^^ 7"!+*, (b > 0). Thus E/T scaling with a non-trivial exponent strongly suggests 
that the underlying physics of the heavy fermion quantum critical point is governed by 
universal physics with du> 3. 



4.4. Failure of the Spin Density Wave picture 



The standard model of the heavy fermion QCP assumes the non-Fermi liquid behavior 
derives from Bragg diffraction of the electrons off a quantum-critical spin density wave 
(QSDW)[^ |4H9|, |9^. The virtual emission of these soft fluctuations. 



e ^ e + spin fluctuation 
generates a retarded interaction 



(146) 



Xq 



Q)2 + ^^2_^ 



(147) 



between the electrons, where xi^l^ (o) is the dynamical spin susceptibility of the collec- 
tive modes. The damping term — /(y/Fg of the magnetic fluctuations is derived from the 

linear density of particle-hole states in the Fermi sea. and T^^ = Tf.^'^ are the 



QCP AFM 




FIGURE 23. Quantum spin density wave scenario, where the Fermi surface "folds" along lines 
separated by the magnetic Q vector, pinching off into two sepai^ate Fermi surface sheets. 



inverse spin correlation length and correlation times respectively. In real space, 

V^ff{r,co = 0)oc——e''i-' (148) 

is a "modulated " Yukawa potential whose range ^ ~ {P — Pc)^^ ^ oo at the QCP. 
Unlike a ferromagnetic QCP, the modulated potential only affects electron quasiparticles 
along "hot lines" on the Fermi surface, that are separated by the wave-vector Q and 
satisfy e^^ = Cjj^q. At a finite temperature, electrons within a momentum range ~ 



are affected by this critical scattering (Fig. ^.). This limits the ability of this singular 



potential to generate non Fermi liquid behavior. There are then two major difficulties 
with the (three dimensional) QSDW scenario for the heavy fermion QCP: 

1. No breakdown of the Fermi liquid Away from the hot lines, the Fermi surface 
and Landau quasiparticles remain intact at the QCP. Thus the specific heat and 
typical quasiparticle mass do not diverge but exhibit a weaker singularity, Cy/T = 



Yo-AVT in the QSDW picture[|9|]. 
2. No E/T scaling The quantum critical behavior predicted by this model has been 
extensively studied p6|l . In the interaction V^j:j:{q,(o) the momentum depen- 
dence enters with twice the power of the frequency, so 

(z=2). 

In the renormalization group (RG) treatment time counts as z space dimensions 
so the effective dimensionality is Drr = d+z = d + 2. The upper critical dimension 
is set by D^^^ = 4, or du = 2[|96|], so 3D quantum spin fluctuations will not lead to 
E/T scaling. In three dimensions, QSDW theory predicts that the scale entering 
into the energy dependent response functions should scale as T-^/^, with a non- 
universal prefactor[| 



4.5. Towards a new understanding. 



In this last lecture, I would like to give you a sense of the seriousness of the failure 
of the spin density wave scenario and share with you some of the new ideas that are 
circulating. Some have argued that it may be possible to explain the of E/T scaling 
and the logarithmically divergent specific heat[|97|] by supposing that the spin fluctua- 
tions form a quasi-two-dimensional spin fluid [pl[ pT] ], lying at the critical dimension. 
Inelastic neutron scattering experiments on CeCu^ ^u^, (x=0.1) support a kind of re- 
duced dimensionality in which the critical scattering is concentrated along linear, rather 
than at point-like regions in reciprocal space[|7^, |9^]. More recent data[|98|] may support 
quasi-2D spin fluctuations at intermediate scales in CeGe2Ni2- 

The assumption that the spin fluid is two dimensional is hard to reconcile with the fact 
that the developing order is fully three dimensional, and with the fact that these systems 
exhibit very little dimensional anisotropy. Even if we accept these problems, there other 
difficulties. First- quasi-two dimensionality can furnish E/T scaling, but it does not 
drop the theory below its critical dimension, and hence has no way of accounting for the 
anomalous exponents in the E/T scaling. 

Finally, there is another more serious difficulty. It has recently become possible to 
examine the approach to the heavy electron quantum critical point through the use of 
field tuning ^ The material YbRh2Si2_j^Gex with x = OA lies precisely at a 
quantum critical point. By applying a small magnetic field, this system is driven back 
into the Fermi liquid. As the field is reduced and the system is tuned back towards the 
quantum critical point, the A coefficient of the resistivity is observed to diverge as 

A oc — 

B 

Such behavior can be obtained in a two dimensional spin fluid model in which the inverse 
squared correlation length is assumed to be proportional to B, ^ ^ oc B. The same model 
predicts a weak dependence of the linear specific heat on magnetic field 

%h°^Logil/B) 

so that the ratio 

Tf, BLog{\/By 

The same experiments also show that the linear specific heat diverges much more rapidly 
with 5, as 70c so that the Kadowaki Woods ratio 

A/'f- = constant. (149) 

It is difficult to understate the importance of this new result. The constancy of the 
Kadowaki Woods ratio over more than a decade in 7 indicates that the momentum 
dependence of the scattering amplitudes in the Fermi liquid are not radically affected by 
the magnetic field, as they would be if the chief mechanism for the mass renormalization 
were derived from the exchange of soft magnetic fluctuations in a 2D spin fluid. These 



new results can only be understood if, in the approach to the quantum critical point, 
the Fermi liquid scattering amplitudes remain local, depending only on the size of the 
renormalized Fermi temperature r^[^. 




E 

ip 



(150) 



In the previous chapter, I argued that the effective Fermi temperature of the Kondo lattice 
measures the "phase stiffness" associated with the amplitude to form of a composite 
heavy electron, so that 

The constancy of the Kadowaki Woods ratio, the lock-step divergence of both A and 
7^, and the appearance of local features in the spin correlations at the quantum crit- 
ical point are all in keeping with the idea that the Kondo bound-state phase stiffness 
is going to zero on the paramagnetic side of the heavy electron quantum critical point, 
just as the spin-wave stiffness goes to zero on the magnetic side of the same point. 
In other words 



> quantum bicriticality ? (151) 



In other words, the Kondo composite bound-state appears to die at exactly the same time 
magnetic order develops. This strongly suggests to me, that perhaps the heavy electron 
quantum critical point might be better understood as a quantum bicritical point, where 
two order parameters go to zero at a point. 

Traditionally, theories of phase transition are built apon an underlying mean-field 
theory. The spin density wave scenario is a consequence of examining fluctuations about 
the Stoner and Slater view of itinerant magnetism. If this approach fails, then perhaps it is 
a sign that we should search for a new kind of mean-field theory to describe the quantum 
phase transition between antiferromagnetism and the heavy electron fluid. There are two 
kinds of suggestion that have been considered recently : 

• Local quantum criticality. The apparent momentum independence of the localized 
critical correlations at the quantum critical point [[7^ has led to the suggestion that 
the the correct mean-field theory, is one that is local, yet fully dynamical. QlOOj 
102[ ] Such "dynamical mean-field theories", [ |103[ ] are thought to asymptotically 



exact in infinite dimensions. In this philosophy, the local physics remains strongly 
interacting even in infinite dimensions, but the local character of the interactions is 
supposed to be stable against finite dimensionality. This idea forms the basis of a 
recent theory by Si et al. 

Traditional RG approach on a new Lagrangian. Rather than abandon the traditional 
RG approach first suggested by Hertz, we should continue to embrace the 
notion that a Wilsonian approach, where interactions become weak in high enough 
dimensions does work for quantum critical points. This approach argues that what 



is needed, is a new description of magnetism, and the way it couples to the Fermi 
sea. One idea here, is that at the quantum critical point, the heavy electron breaks 
up into its spin and charge components. [|68|] 

We now discuss these ideas in more depth. 



4.6. Local Quantum Criticality 



The momentum-independent scaling term in the inverse dynamic susceptibility (7) 
suggests that the critical behavior associated with the heavy fermion QCP contains some 
kind of local critical excitation [[72||. One possibility, is that this critical excitation is the 
spin itself, which would then presumably develop a slow power-law decay [ ]KR]| , |1U1 



102] 



(5(t)5(tO) 



1 



(152) 



where e 7^ signals non- Fermi liquid behavior. 

Si[ |104| | et al. have extensively developed this idea, proposing that the local spin 
susceptibility Xioc = ^)la)=o diverges at a heavy fermion QCP. From ( |143| ), 



d'^q- 



1 



q_Q)2 + jc 



j(d-2)a/2 



(153) 



so a divergent local spin susceptibilty requires a spin fluid with J < 2. Si et al are thus 
motivated to propose that the non-trivial physics of the heavy fermion QCP is driven by 
the formation of a two-dimensional spin fluid. Si et al consider an impurity spin within 
an effective medium in which the local Weiss field H has a critical power-spectrum (Fig. 
H) 

(|//(a))|2)^;to"i((o) = «^ (154) 




e-l 



X-\u) = {\H{u)\') 



X 



-1 = Xo' + ^4a;'''^ 



FIGURE 24. In the local quantum critical theory, each spin behaves as a local moment in a fluctuating 

Weiss field. In the theory of Si et al [102], a self-consistent solution can be obtained for e = 1 in which 

1 

the local susceptibility develops a self-energy with a non-universal exponent. M(a)) oc opx. 



where e is self-consistently evaluated using a dynamical mean-field theory, where q— 
dependence of self-energies is dropped. In principle, the method solves the dynamical 
spin susceptibility of the impurity X^^{^) = Xo^{(o) +M{(o). This, in turn furnishes 
a "spin self-energy" M{(0) used to determine the spin susceptibility of the medium 
X~\q,co) =J{q)+M{C0). 

Si et al find that a self-consistent solution is obtained for £ = 1, if the spin-self energy 
contains a separate power-law dependence M{(o) ~ ft)" with an exponent a = 1/pA 
which is determined by the density of states p and band- width A of the bond- strengths in 
the two-dimensional spin fluid. Although self-consistency requires a new power-law in 
the spin susceptibility, independent solutions of the impurity model have not yet shown 
that this feature is indeed generated by a critical Weiss field. This theory nevertheless 
raises many interesting questions: 

1. Is the requirement of a two dimensional spin fluid consistent with the ultimate 
emergence of three dimensional magnetic order. For example- does the the cubic 
(and hence manifestly three dimensional) quantum critical material, Celn-^ display 
a divergent specific heat? 

2. If the spin-fluids are quasi-two dimensional, do we expect an ultimate cross-over to 
a three-dimensional QSDW scenario? 

3. If a is non-universal, why are the critical exponents in CeCu^^Aux and 
YbRh2Si2_^Gex so similar? 

4. What stabilizes the local quantum criticality against intersite couplings? 

4.7. Ideas of spin charge separation and supersymmetry. 

An alternative possibility, is that the heavy fermion QCP is a truly three-dimensional 
phenomenon. In this case a different approach is needed- we need to search for a 
new class of critical Lagrangian with du > 3[ |1U6| ]. On general grounds, the existence 
of a Fermi liquid in the paramagnetic phase suggests that the new class of critical 
Lagrangians must find expression in terms of the quasiparticle fields i// in the Fermi 
liquid- but how do we couple these degrees of freedom to the magnetism, and how do 
we account for the simultaneous loss of the resonant bound-state stiffness at the same 
time that magnetism develops? The simplest possibility is to write 

L = L^[xif]+ L^_j^[xff, M] + [M] . (155) 

where Lp describes the heavy Fermi liquid, far from the magnetic instability, de- 
scribes the magnetic excitations that emerge above the energy scale (P) . 

Lp describes the way that the quasiparticles couple to and decay into critical 
magnetic modes; it also determines the type of transformation which takes place in the 
Fermi surface which occurs at the QCP. This last point follows because away from the 
QCP, magnetic fluctuations can be ignored in the ground-state, so that Lj^ 0. In the 
paramagnetic phase, (M) = so — > 0, but in the antiferromagnetic phase (M) ^ 0, 



I.e. 

^f[V^] paramagnet 



^eff 



L}[¥]+hM[wAM)] a.f.m. 



where the asterisk denotes the finite renormalizations derived from zero-point fluctua- 
tions in the magnetization. 

If the staggered magnetization is the fundamental critical field, then we are forced to 
couple the magnetic modes directly to the spin density of the Fermi liquid 

k,q 

But once the staggered magnetization condenses, this leads directly back to a static spin 
density wave (Fig p3t ). 

An alternative possibility is suggested by the observation that the magnetism develops 



spinorial character in the heavy Fermi liquid. The Luttinger sum rule[ 108] governing the 



Fermi surface volume V^^ " counts" both the electron density and the number of local 



moments per unit cell n^^^^^^]^, ^ 



2(1^3 ="^ + %».- (157) 

The appearance of the spin density in the Luttinger sum rule reflects the composite 
nature of the heavy quasiparticles, formed from bound-states between local moments 
and high energy electron states. Suppose the spinorial character of the magnetic degrees 
of freedom seen in the paramagnet also manifests itself in the decay modes of the heavy 
quasiparticles. This would imply that at the QCP, the staggered magnetization factorizes 
into a spinorial degree of freedom M{x) = z^{x)dz{x), where z is a two-component spin 
1/2 Bose field. "Spinorial magnetism" affords a direct coupling between the magnetic 
spinor z and the heavy electron quasi-particles via an inner product, over the spin indices 

L?-M = SBX^^ (^^-qaV^ka) +H.C], (158) 
k,q 

where conservation of exchange statistics obliges us to introduce a spinless charge e 
fermion X- This would imply that the composite heavy electron decays into a neutral 
"spinon"and a spinless charge e fermion ^ S(y ■ The critical Lagrangian in this 
case would take the form 



L = L^[v/] +Lp_^[\{r,x,z] yr]. (159) 

We have to be cautious of course, because this is undoubtedly one of many alternative 
ways we might begin to construct a new class of critical Lagrangians. What we do see 
quite clearly however, is this line of reasoning leads us to into the notion that the break- 
up of the heavy fermion QCP involves spin-charge separation. 

Hall constant measurements may provide a good way to discern between the spin 
density wave and composite quasiparticle alternatives. In the former, regions around the 



hot-line do not contribute to the Hall conductivity, and the change in the Hall constant 
is expected to evolve as the staggered magnetization . By contrast, the composite 
fermion scenario leads to a much more rapid evolution: provided that the density of 
spinless fermions is finite at the QCP the Hall constant will jump suddenly at the QCP 

I Mr,, (vectorial ) 

^«°=|o(T) (spincrial) ("^O) 

The only available Hall measurement at a QCP to date shows a change in sign takes place 
in the close vicinity of the QCP in critical CeCu^ ^u^, it is not yet clear whether there 



is a discontinuity at the transition 01O9| ]. This is clearly an area where more experimental 
input is highly desirable. 

Let me end with a few speculations. If we are to construct a new critical theory for 
the heavy electron quantum critical point, then we will need some new theoretical ideas. 
A new critical theory will require as a first step, a new kind of mean-field description 
that permits us to understand why the magnetism and the Kondo effect die at a common 
critical point. At present, we do not know how to construct a mean-field theory that 
contains a heavy electron quantum critical point. One interesting idea here, may be the 
incorporation of supersymmetry. This is an idea that I have tried to develop, as yet 
with only partial success, with my graduate student John Hopkinson and collaborators, 
Catherine Pepin and Alexei Tsvelik.[[106|, |107|] At a very naive level, magnetism involves 



the manifestation of spin as a bosonic excitation, whereas heavy electron behavior 
involves the manifestation of spin as a fermionic object. If the two phenomena share 
the same quantum critical point, then is it possible that the spin manifests both types 
of behavior at a quantum critical point- in otherwords, that it displays some kind of 
supersymmetry? 

It does prove possible to represent both spin and Hubbard operators in a way that 
involves a locally supersymmetric gauge theory, but a mean-field theory still eludes us 
at the current time. We do have an idea about the structure of this mean-field theory, 
which I shall briefly mention to you. In this putative mean-field theory, we require two 
order parameters- one for the formation of the composite heavy electron corresponding 
to the amplitude for composite fermion formation ~ V and one for the magnetism 
(V/2 ~ (^tr) ~ VM). Suppose we can integrate out all of the fermions in the theory, so 
that we are left with an effective theory for and i/Zj, given by a Landau Ginzburg 
free energy , i//2] • Now here's the remarkable thing- for the two order parameters to 
share a quantum critical point, then the expansion of the Free energy near the QCP must 
take the form 

F~a(|i//2|^— |V^il^) + interactions 

When a > 0, we have the heavy electron phase with i//2 = 0, 1//^ 7^ but when 
oc < 0, we have the magnetic phase, where 7^ 0, = 0. At a = 0, both order 

parameters vanish simultaneously. The two must go to zero at the same point in the phase 
diagram, so they can only come together in the quadratic combination (Iv^jP — |V^2l^)- 
This suggests that the negative definite metric 

(IV^2l^-|V^iP) 



Spins ~ bosons Spins ~ fermions 




FIGURE 25. The heavy fermion QCP may involve a supersymmetric gauge symmetry. To understand 
the fact that the magnetization (M) disappears at precisely the same point in the phase diagram where the 
amplitude (V) for the formation of composite fermions goes to zero, we need a special symmetry between 
the two order parameters. Local moments behave as fermions in the paramagnetic phase, but as bosons in 
the antiferromagnet. Does a supersymmetry develop at the quantum critical point and is this responsible 
for the "Minkowskii" metric between the two order parameters that is required for them to vanish at the 
same point in the phase diagram? 

is a symmetry invariant of the quantum critical point. The appearance of a minus sign 
- a Minkowskii type metric- is required by the phenomenology, yet traditional invariant 
symmetry groups of a critical point involve a positive metric associated with a trace over 
order parameter combinations. The minus sign would occur if the residual symmetry 
between these two order parameters corresponded to symplectic group. One way for 
such minus signs might appear is via the supertrace- an invariant of a supergroup. This 
prompts the following conjecture, on which I will end: that the the above invariance is a 
residue of a supertrace in a supersymmetric Lagrangian for quantum criticality. 

4.8. Summary 

This section has discussed the origin of the mass divergence at the heavy fermion 
quantum critical point, emphasizing that a quantum spin density wave picture can not 
explain the observed properties. The proposal of fundamentally new kinds of quantum 
critical points has been reviewed. This is clearly an area with a huge potential for 
progress both on the experimental, and theoretical front. 



4.9. Exercises 



1. Consider the tree level scaling for the zero temperature Hertz-Millis Lagrangian at an 
antiferromagnetic quantum critical point 



S = S^ + Si (161) 



where 



IcoPA 



{k-Qy + 



(162) 



describes the propagation of an overdamped spin fluctuation with inverse susceptibility 
X~^(k,(o) at quantum criticaUty, with a critical wavevector Q. Here we have used the 
notation K = {k,(o) to denote the wavevector and frequency of a magnetization mode. 
The non-Unear interaction term in the model takes the form 

^iiu]=u J ( .n ^2^)^^''i^---^^''4^^'^'^£''^-^ ^^^^^ 

(a) Derive the tree level scaUng that keeps the Kinetic term 5^ invariant. First note that 
from the kinetic term, the scaUng dimension of frequency is [co] = [kY. Show that if 
the wavevector and frequency cut-off are rescaled according to 



A ^ CO, 

A = 

then one must rescale 



A=--, cd, = -^ (164) 



k = kb, 03 = &)V (165) 
to keep the form of the Kinetic term invariant. Show that under this scaling 

^ Aa)\<a„,k<h ^d+lj^ 
^ b(^+1-2) j -^^iX-' ( mi K)M{-K) 

SO that with scaUng 

M(jc)=M(k:)&-(^+''-2)/2 

the kinetic energy remains invariant, 
(b) Using the tree level scaling derived above, show that the interaction term transforms 
as 

Si{U)^SriU*) 

where 



^,(z+<i-4) 

showing that the interaction term scales to zero for <i + z > 4, proving that d = 4 — z- 
2 is the upper-critical dimension for the Hertz-Millis model. 
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